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SECTION  1 

AMBIENT  ATMOSPHERE  MODEL 


INTRODUCTION 

The  ambient  atmosphere  models  used  in  the  current  WEPH  code 
are  essentially  the  same  as  those  developed  by  Science  Applications 
Incorporated  (SAI)  for  ROSCOE  (References  1-1  and  1-2).  Quantities 
modeled  include  pressure,  density,  temperature,  density  scale  height, 
and  major  and  minor  neutral  species. 

The  altitude  region  is  divided  into  a low  altitude  region 
(h  < 120  km)  where  the  major  and  inert  species  are  thoroughly  mixed 
so  that  the  fractional  concentrations  arc  (almost)  altitude  indepen- 
dent, and  a high-altitude  region  (h  > 120  km)  where  diffusive  separa- 
tion prevails.  In  the  low-altitude  region  a day-night  change  in 
certain  species  concentrations  is  provided  for  but  there  are  no 
changes  with  season,  latitude,  or  solar  cycle.  In  the  high-altitude 
region  the  species  and  properties  are  modeled  as  a function  of  solar 
flux. 

Recently,  SAI  revised  the  ROSCOE  atmosphere  model  for  use  in 
optical  predictions  (Reference  1-3).  The  revision  includes: 

1.  Addition  of  new  species;  N(2D),  N^O,  CO,  CH4,  O^D), 

H,  OH,  H02 

2.  Provision  for  determining  atmospheric  temperature  be- 
low 120  km  as  a function  of  latitude  and  season 

3.  Provision  for  determining  the  species  0,,  H-O,  N-,0, 

4 2 ill 

N(  S),  N(  D),  and  NO  as  a function  of  latitude  and 

4 

season.  Provision  for  determining  the  species  N(  S) , 

2 

N(  D),  and  NO  as  a function  of  local  time. 
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Most  of  the  above  changes  have  been  incorporated  into  a re- 
vised WEPH  code  atmosphere  model  (the  species  CO  and  CH^  are  not 
modeled).  In  addition,  the  structure  of  the  revised  WEPH  code  atmos- 
phere model  and  the  smoothing  procedures  used  were  changed  to  facili- 
tate use  in  propagation  codes. 

In  the  WEPH  Code  atmosphere  model  (and  also  in  the  ROSCOE 
model)  initialization  calculations  arc  made  for  a given  location, 
date,  and  time  prior  to  determining  atmospheric  properties  at  speci- 
fic altitudes.  The  initialization  calculations  require  use  of  a 
number  of  special  routines  that  arc  only  used  during  initialization. 

In  the  current  model  the  initialization  calculations  and  the  atmos- 
pheric evaluation  calculations  are  made  from  the  same  routine  result- 
ing in  all  the  atmosphere  routines  being  present  (in  core)  when  atmos- 
pheric properties  are  calculated.  In  the  revised  atmospheric  model 
the  initialization  calculations  and  the  evaluation  calculations  are 
separated.  This  allows  overlaying  the  initialization  routine  to 
reduce  core  storage  when  initialization  is  only  done  once  for  a given 
region. 

A major  objective  in  the  development  of  the  atmosphere  model 
by  SAI  was  obtaining  continuous  derivatives  for  modeled  quantities. 

To  this  end  polynomial  fits  are  used  to  describe  quantities  between 
specified  altitudes.  The  use  of  high  order  polynomials  (12th  degree 
in  many  cases)  can  lead  to  numerical  difficulties  in  solving  for  the 
coefficients  and  in  subsequent  evaluation  of  the  polynomial  for  some 
machines.  As  described  in  Section  2,  Booker  (Reference  1-4)  has 
suggested  use  of  an  exponential  formulation  for  fitting  multi-region 
ionospheric  profiles  of  electron  density  that  provides  continuous 
derivatives.  The  formulation  suggested  by  Booker  (sec  Equation  2-10 
in  Section  2)  can  be  used  as  an  interpolation  procedure  by  choosing 
the  smoothing  scales  so  that  the  effects  of  smoothing  outside  the 
interpolation  region  are  negligible  within  the  interpolation  region. 
For  the  atmosphere  model  the  smoothing  scale  was  chosen  as  one-tenth 
the  distance  between  data  points. 
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For  this  choice  of  a 4-point  interpolation  can  be  used.  Let  N(z) 
be  specified  at  Zy  z^,  Zy  and  z^.  Then  N(z)  between  z?  and  z^  is 
given  by 


£n[N(z)]  = fn[Nr]  + A^Zj  " zy) 

+ (A2  - At)  j f ( z - z2,  B2)  - f(zy  - z2,  B2)j 

* (A,  - A2)  }f(z  - zs.  B,)  - f(zr  - z3>  Bj){  (1-2) 


Figure  1-1  shows  values  obtained  by  interpolation  between  data  given 
at  5 km  altitude  intervals.  Note  that  the  interpolated  values  do  not 
exactly  reproduce  the  data  points  since  this  is  where  the  smoothing 
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INTERPOLATION 


Figure  1-1.  Illustration  of  profiles  obtained  from  interpolation  model 


is  done  to  provide  continuous  derivatives.  If  the  derivative  of 
N(z)  is  required,  it  can  be  obtained  from 


dN(z) 

dz 


N(z) 


I A A2  ' A1 

| A1  + 1 + exp[-  (z-z2)B] 


d-3) 


A3  ' A2  ) 

1 + exp[-  (z- 7 B]  j 

In  the  WEPH  code  (and  ROSCOE)  chemistry  models  effective  photo- 
dissociation rates  for  photodissociation  of  C>2,  0^,  and  NO^  are  deter- 
mined from  ambient  species  concentrations.  As  discussed  in  Section  4 
the  use  of  the  nominal  SAI  model  for  NC>2  results  in  a negative  effec- 
tive photodissociation  rate  for  N02  above  about  50  km.  In  order  to 

prevent  this  the  ambient  daytime  NO*  is  limited  to  a value  determined 

^ -3-1 

from  a nominal  photodissociation  rate  (6.8  x 10  sec  ) used  by 
Schicbe  in  the  DCHEM  code. 

Tables  1-2  and  1-3  show  results  obtained  with  the  revised 
atmospheric  model  for  conditions  shown  in  Table  1-1  (see  Section  6 
in  Reference  1-1  for  values  obtained  with  the  current  atmospheric  model). 
Potentially  significant  changes  to  species  concentrations  previously 
modeled  are  the  reduction  in  NO  and  H20  in  the  D-region.  The  NO  con- 
centration has  been  reduced  by  about  a factor  of  5 and  the  H20  con- 
centration has  been  reduced  by  a factor  of  2 in  the  lower  D-region 
to  about  a factor  of  5 in  the  upper  D-region.  Results  from  the  ambi- 
ent ionosphere  model  are  also  shown  in  Tables  1-2  and  1-3.  The  quan- 
tities VEM,  VIM,  and  VEI  shown  in  Tables  1-2  and  1-3  are  the  electron- 
neutral,  ion-neutral,  and  electron-ion  collision  frequencies,  respec- 
tively. Changes  to  the  ambient  ionosphere  and  ion-neutral  collision 
frequency  models  are  described  in  Section  2. 
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Table  1-1.  Conditions  for  Tables  1-2  and  1-3. 


Table  1-2 

Table  1-3 

North  Latitude  (deg) 

45 

45 

East  Longitude  (deg) 

235 

235 

Year 

1977 

1977 

Month 

September 

September 

Day 

1 

2 

Hour 

12:00 

0:00 

Solar  Flux  (x  1022  W m"2  Hz-1) 

158 

158 
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Table  1-2.  Revised  ambient  daytime  atmosphere 
species  and  temperatures. 
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Table  1-2  (Continued) 
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Table  1-3.  Revised  ambient  nighttime  atmosphere 
species  and  temperatures. 


f 


I 

\ 


ov  w 

2U  UiUi  tA 

x »>5  Ox 


X 

ml 

is 


is 


3B 

W x 


U 

* 


-oo2o82o-2o22o22o222--o--o2-oo-oo2oS2o2oo2SoSooS 

11IIII9IIIIII  911 

wiiiMMkiwmiitHi1  wwwwMtf  WWWWVIW 

•<**K«ehf«Q|«te«NO»Ne#NO«MonriON<«ONOo*«ooNoooooofo«aeN 

— 00-60  — oo  — ©o  — oo-oo- 00— 00-00-00  — 00  — o c — 00  — ooccoooo 

^ ^ ^ yi  ^ ^ U|  ^ ^ ^ ^ yj  y|  ^ yj  yj  yj  yj  yj  yj  yj  y^  yj  yj  ^ y|  yj  y^  yj  yi  y|  y|  yj  y|  yj  y|  y^  yj 

— o — «4»m««»«  — — o?a-*«+M©i‘MMO^in««MM»««»»'fttf>M'rM 

»Omamo*«om^4Nm«4om^iraAiN«^x»«omKm«««aoaN^MO»h«omavt 

-So-Sc25S22o2SSooSooSSoSooocSo£oS£5c2£ocSooS*SSS 

iiHimMiiimilnimiMiMiminimMimi 

y^  ^ yj  y^  y|  y^  ^ ^ yj  y^  y^  ^ ^ yj  ^ y^  ^ yj  yj  yj  yj  y^  y^  y^  yj  ^ ^ y|  yj  yj  yj  ^ ^ yj  y^  yj  u yj  ^ ^ yj  yj  y^  |y  yj  yj  y^ 

«m«  — mm*o««oo«m«  + — «o»’"in»M«M<**— «r«u— m»o  — «*  — ^ftO**©***© 


■-  n 1 


i n 

“ 4U 
< SJ 

ou  ................ 

^ ^ y|  yj  yj  yj  yj  yj  yj  yj  yj  ^ yj  yj  yj  yj  ^ yj  ^ yj  yj  yj  yj  yj  yj  y^  yj  yj  yj  yj  yj  yj  ^ yj  yj  yj  yj  yj  yj  |y  yj  ^ yj  yj  yj  yj  yj  yj  yj 

ll  3s2s5si22js25?2552s:s2isssssssissss;ss5;ssps:s5rs 


X w msj  <j 

QV  IV  U 

? X X Ml  X 


_ u ou 
OU  X u u 

t x X OU 


oo^Sooo^oo^oooo^oo^ooSooooooooooo^ooSo^ooooo 


i*  2 


m • 

0 

0 

§ 

g 

U «N  • 

O 

C 

• • 4 N 

0 

0 

h © — * 

O 

O 

«4N* 

g 

« 4ii  i 

O 

O 

pt  00 

•© 

>■ 

itfs  M - 

0 

C 

• m«i 

0 

0 

l©V< 

© 

© 

i 4 •+* 

O 

O 

• ♦m- 

0 

© 

* 4 f*  4 

O 

O 

> tr  — f 

g 

Ui 

Ui 

Ui 

3 

O 

• 

— 

Ui 

Ui 

Ui 

Ui 

Ui 

Ui 

Ui 

Ui 

i.ooc 

8 

• 

— 

§ 

• 

g 

• 

g 

• 

g 

• 

O 

O 

• 

c 

0 

• 

O 

O 

• 

g 

• 

g 

• 

— 

g 

• 

g 

• 

g 

• 

O 

O 

• 

— 

* 

m «r 

1**  r> 

<n  « 

m 4 

m 4 

**  t* 

m ** 

m m 

«M 

m m 

m m 

« M 

m — 

t+  — 

— 

0 — 

0 — 

0 — 

O — 

O — 

0 — 

0 - 

0 - 

C - 

c — 

O — 

© — 

© — 

© — 

0 - 

Ui 

u ui 

Ui  Ui 

Ui  Ui 

Ui  Ui 

*¥ 

Ui  Ui 

UiUi 

UiUi 

KS 

Ui  Ui 

Ui  Ui 

'it 

St* 

Ui  Ui 

** 

IT 

^ 4 

***  M 

I**  © 

«r  ir 

**  K 

il 

c - 

0 - 

C M 

• « 

+ • 

© e- 

► ** 

M — 

M ^ 

1*  M 

« • 

► © 

• *• 

© c 

© M 

M — 

4-  M 

tr  — 

• 

* ♦ 

— M 

#*  — 

IT  « 

« M 

© — 

© * 

♦ M 

m — 

M* 

M 4 

— M 

— — 

- « 

**»  m 

om«»om»»»m«£— *r-<r?o  + £C+£»  + »««**^*©p»m««M«v«Mt»*  — 

--0--00  — 00-00-00— oo-ooocoeoooooooooooooooooocooo 

itiiiTiiTiiiii 

yjyjyWWylylyjyjyyjyjyjyjyjyjyj^yiyjyjUyjyyjy,y,y,yjy,yjyjyjyjy|yjyjUyyjyjWyjWy'yjy|yj 

o--cNo»#»H<»>^M«HN40NfON4c«'iro»^-wir-«r«^Nrrx^«'*r^c 

ooocoooooooocoocoooooccocecccoooeoooocooocoooocc 
I • • I I • I 

yyuyyyyyiyyyyyyyyyyyyiyyyyyyy 

f*  O m *■  «»0<A*M*>»0P»*»«M  — O <*  « * m « — * 


I 4 O « 4 • « — « — • — 


U XU  — u 

U OW  u*  Si 

C X X > x 


U U W W 

MW  XU 

OX  X > X 


C u 

U — U lw 

MU  — U W Si 

* X MX  > X 

— o — — 


eMO0M»ONheM^CM40N«0««0ff««0fll«0«f)0«ir0«#Cli0#lKp«N««*« 

co  — cooooooooooooooooocooooeooocoooeoooooooooooec 
• I • I • • I • I • I • • • • I 

y y y y yi  y|  yj  y yj  yj  yj  yj  y yj  y y y y y y yyyyyyyy  yj  y y y y y y y y/  y ^ y yj  y y y y y y y 

o o — c — «o  — ‘*o«*o«ec*«ot»MO*Mc^r<e«»c»cc<"  + o«  — o«nM^«^*M  — 

-►*-•*  — — O — *■  * — M f-  — MM  — — — — M — — ►»—  M * — « © ► «—«■ 

_ M — - «.*M  — — M*—  « — — «*  — N©— • — — *M  — M*  — — O — «— — «*  — *M«  + ♦ — *»  ► 

• o»*o»«otho«^c«hOt«^acK«o^iro«<rc«tf«cfiro««Q««cir«Nf 
-00— 00  — 00  — 00  — 00  — oo  — oc-oc-co— oc  — 00  — co  — 00-00— 00  — 00 

y y y yj  y y y y y y yj  y y y y y y y y y y y y y y y y y y y y y y yj  y y y y y y y y y y y y y y 

K C O *■  O S"  O O O #•  C <N  .T  C * C O X C « ^ G * - CKNiW^C 

*00— 0*©C«©O«C0r«O©«0*©O»«0»M0©©0M40M«00©0M— iTOC*© 

^o  — »oo«oo«oo«oje*or^c^^c««oc*«-««oo^c^irc^^©e“rf'c#^M< 
— o — — o — — o — — o — — 00  — 00  — c 0 — 00  — 00  — CO  — 00  — OU  — oc  — 00  — o o — c o 

Uttiuuwit/iMiuUwyuyyyyyyyyyyywywywwyyyyyyywyyyyyyyyyyw 
»0<*©0M*O  — ©O-OOMOCM*-  CO*0**OMOO—  • 0»»04*>0*00«<V*’U«ON 
► 4 * - # M 4 4 — M + » * « p»  « + o - 4 »«#—•*««  * * f>  + M t*  — 


§ s 


s 


o o 

• • 

o ^ 


S i 

• • 

* o 


19 


I 


ou  io  i *t 

♦V  O W M 

X V.  *-  o o >. 


ft  — — ftft  — ftN  — *•►■**  — 


’«  * 


► u 
►*  o 

— »*\  M O 

O U *»  «A  < M 
MU  * * J •/> 

— x i*. « a v 

*-  °5 


r««4»affi«»tAOi 

ooooooooocooo  — o 

< I I • I I I « I I I • I I ( 

<v<?  — 

• — o 400  airoomoM^c 


• «coft«of^ft  — *>  « — *• 

>0000000000000 
• 1 1 11  < 1 


fft  — m-  o — 

-*oft©*ocoi*fto  — ft 


« u * x X 

— x,  w x a n 

X — ft  ► -J  o 

So  <0 


NU  O 

CO  ♦ o 

O V f \ 


WJ  uw  u>w  u w w u 

♦ *•  o «r  — « ci 

♦ ac  ft*  »*  ft  «« 

• • • t I • I • 

— -ft  -IV  M • N » 


MU  *Nt  O O 

SO  X O O 

•s.  ^ W v 


>0000000004 

JWWWWWWWWWi 
MV-COft- OOC  I 


'44C4#CNMC« 

;ooooocooo< 


> IV  » * ft  ft  *»  *>  -O  I 
14  ftWMr^-rf'ar 


o X O — w 

O CO  U.  4) 

C N V > >. 


u o t w 

IV  o X o — u> 
O'*  N > S 


—00—00—00—00—00 

W MWMMWWWMWMWMMM 
— CftCO  — 4»mMftl>NC4' 
irfvftftftfvrf'M  — f»  ft  « — Oft 


IV  O — O w/ 

X "V  IV  ^ > "V 

— 0 — — 


20 


Table  1-3  (Continued) 


*8 


i £8 


28 


oin  — ANgAfNOr»NQANQi»N  — an  — AN  — AN  — nn— A'NNANNONNONiNOiNr*or«m«A* 

O OOOOdOOSo  0 60  06  0 0^0  0^0  0(j0  0^0000  0(|0  0^6  00d0lj0  0(j0  0(j6  0<j 

“^*!J^**SS*y«»gtfrSSS!S!Sf!K^lS5!tgSSSf5iSJiiS»aiS|!t»5;SSSSSSS?P^«« 


z 

% 


a sj 

s* 


i! 


H6oH5^6o(jiO(j6oooo^6oijOoo6eo5oooooooooo(joooooo6coooo 

aa*a**S8**8***sssa?«a5[«^«a*sas**»tt»»*B»asa*«ssaa*a* 

-eNOAAeAOOOANOA  — — aoa  — AAOeAAooeo-*— oooaaaoao  — noainaon 


— N-  --  ON  ON  ON  ON  ON  OA  ON  ON  O A ONI  ON 

— Q - ON  O - O - On  On  On  On  O n 6 n On  On  5 — O - 


17  17  17  17  17  17 


4 


* V 

* U 

— *s  ... 

« - ui  * 


28  88  88  88  88  88  88  88  88  88  2$  88  88  88  88  88  88 

' ••  • 

— NN 

i:i$i 


oooo^ooooo^oo^oo^loo^oo^oS^oSooSSso^oc^oo^ooooo 


is 


o u 

— V 

— N» 

o — 


N u 
ou 
u -v 


SJ 
N U 
X >* 


On 

i/CC 

8888888! 
N # # N ► O O 

iiiioN 

JU  Ui  i 

222222282222' 
b >4  N N**N  NmNNnN 

i$s$?£issi;ss8$sH 

ag 

O A 
• • 

O 

o 

NN  N|  * 

1* 

- a m 

Si 

♦ NN  < 

si 

« 

N -*N* 

1* 

- — O' 

1* 

a a a 

f* 

>»  «i  «i 

1* 

S8  SS  38 

A A A 1 

I8 

N A o • 

SS 

-N  « « 

SS 

-N  A • 

S2 

S8  38 

AO 

A 

e 

SJ 

8 

• 

A 

Si 
• • 

si  3* 

• • k • 

• — IN  — 

si  *i 
• • • • 
♦ i Ni 

= i $g 
• • • • 
e — a a 

ri  5i  si  H si  *!  si  *g  IS 

• >4  A A llix  N*«  oirt  4 — •«••  N •<  A A 

8 

• 

e 

So 

Si 

INN. 

ee 

INN* 

ee 

IN  N 

e e 

ss 

ss 

Af  A«  AA 

ee  ee  ee 

SS 

ss 

ss 

ss 

A « A A 

o o o o 

i 

U 

u 

8 

88 

88 

88 

U Ui 

Ui  Ui 
N*  • 

88 

88 

88  82  88 

5 A 

Mki 
0 — 

28 

88 

Ui  Ui  Ui  Ui 

U> 

A 

Ui  V 

IN 

• 

ON* 

• • 
- 4 

• O 

• • 

e 4 

• • 
N A 

• a 

• • 
— IN 

A IN 
• • 

NN 
• • 
•i  N* 

• • 

— 0*  AO  A# 

-IN  --  — • 

• • 

A A 
• • 
AN 

A N 
• • 
N A 

A A 
• • 
N A 

4 A « 7 
• • • • 
AN  AO 

e 

• 

A 

z 

a 

IN 

o 

IS 

is 

iS 

8S 

••N# 

— IN 
?° 

Is 

••  oe  oo 

A A 
1° 

?• 

I8 

88  8^ 

O 

u 
o w 

V 

8 

• 

— * 

•S 

• • 
- 4 

n 
• • 
- 4 

UJUi 

N a 

4 m 

• • 

— IN 

88 
— o 
• • 

88 
► ♦ 

• • 
*- 

88 
A A 
• • 
♦ — 

88 

• Ok 

NO 

88  8|  ll 

Si 
• • 
AN 

is 
♦ — 

Si 
• • 

4 0 

UiUi 

o o 

• • 
• 1 

28  88 
AO  A A 

• • • • 
A A AO 

Ui 

0 

A 

• 

A 

O 

IN  ► 

o o 

INN* 

o o 

N A 

ee 

IN  A 
OO 

ss 

SS 

SS 

SS  SS  ss 

SS 

So 

N A 

e e 

ss 

N A N A 
O O O O 

A 

O 

t“ 

IN  o 

*"8 

U 

u. 

A 

• 

N 

IN  Ui 

A ♦ 

•N  ♦ 

• • 
IN  IN 

UiU 

4 o 

• • 

IN  IN 

Ui  Ui 

SS 
• • 

NUN 

« 

INN* 

S* 

♦ e 

♦ « 

is 
• • 
N - 

si 
• • 
A — 

si  «« 

♦ - 4 — 1*  — 

UiUi 

A A 
O A 
• • 
N*  A 

Ss 
• • 

A 0 

88 
A N 
• • 

m a 

88 

AO 

• • 
00 

2*  tt 

4 O AO 
• • • • 
O A ON 

Ui 

8 

• 

O — 

O — 

h O 
o a 

- O 

o — 

82 

82 

82 

82  82  88 

8S 

SS 

ss 

8S 

Sc  00 

O 

— MJ 

Ui  a 
> N. 

Ui 

N 

Nl 

IN 

Ui  Ui 

IN  IN 

— O 
• • 

Ui  «i 

e a 

• • 

88 
o in 
i^i 

Ui  Ui 
A A 
♦ A 
• 

88 
— 4 

0** 

Ui  Hi 
N*  N* 
a a 

ON? 

Ui  Ui 
IN  IN 

• P 

■/in 

Ui  Ui  Ui  Ui  Ui  Ui 

A*  A A —A 

o « o A A • 

Ui  Ui 

2n 

•NA* 

Ui  Ui 

=3 

/o? 

88 
a m 

Ui  Ui 
A N 
♦ A 

A A 

?i  SS 

/r*  * • 

8 

A 

• 

- 

A -• 

o -« 

IN  •> 
O — 

NIC 

O — 

in  e 
o — 

IN  O 
O — 

INO 
O — 

o S 

SS  SS  oS 

oS 

oS 

o 2 

oe 

AO  O « 

0 e 00 

C 

O 

r w 

Ui 

UJM 

Ui  UI 

Ui  Ui 

88 

88 

Ui  Ui 

Ui  Ui 

28 

Ui  Ui 

88 

UiUI 

88  88 

8 

> V 

o 

• 

• 

o a 
• • 

in  a 

O O 
• • 

S3 
• • 
N N 

• • 

• • 
A IN 

« »N 
• • 

IN  IN 

♦ 4 
• • 
IN  — 

A • 

• • 
— o 

• • • • • • 

— 4 A a O IN 

N 4 

• • 
A — 

SS 
• • 
A A 

A 2 
• • 

N « 

O A 

• • 
— 4 

SS  -0 
• • • • 
A N • A 

S 

• 

O 

u 

2.8 
> x 

IN 

Ui 

« 

IN 

• 

A 

4 IN 

O — 

88 
— A 
• • 

»N  — 

♦ a 

o — 

88 

— 4 

• • 

— A 

IN  A 

o — 

Ui  Ui 
Ni  IN 
INO 
• • 

« * 

IN  a 
O N 

Ui  Ui 
N*  A 
• IN 
• • 
a in 

NI  — 
O « 

88 
e ♦ 

• • 
IN  — 

NI  O 
O — 

Ui  t*. 
4 »N 

m m 

• • 
— o 

AO 

o — 

88 
A A 
• • 
•*  4 

AO  INO  INO 

o a o — e — 

4 A N*  A ON 

ON*  IN  A A — 

l • • l • • 

— 4 • A ON 

N O 

Si 
• • 

4 — 

ss 

88 

0 m 
• • 

N O 

SS 

SS 
• • 

N O 

SS 

H 
• • 
— ♦ 

NO  a O 

e e ee 

Ui  Ui  Ui  Ui 

O — O — 

O A A A 

• • • • 

AN  A A 

O 

8 

O 

• 

A 

s 

s 

8 

i 

8 

8 

8 

8 

8 8 8 

S 

8 

8 

8 

8 8 

IN 

O 

o 

A 

A 

o 

N 

A 

* 

IN 

A 

i 

A 

4S 

IN 

ON 

S 

a 

a e e 

0 A * 

AAA 

2 

A 

O 

• 

A 

i 

8 

N 

e e 

N O 

N N 

21 


N - 


Table  1-3  (Continued) 


w u 

MU  So  * 51 

I x UJ  £i  MV 


ooo2ooooooooooooooooooooSoooooooooooooooooooooooooo 

I I • I • « I I I I I I ' I I • I 

Itl  WW  111  til  UJ  U 111  111  U y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y 

-*ix»o<Mmm4om  4 rx  • 4 ir  44m  444  44441x4  44a  444  44  — 4 4 a 4 4 4 4 4 <n  # « mi  « 

404400404400  — 04aC4mo4Mi04ooin404'n04404Mioinoo4a044oa404 


— 4 4 — 4 4 — 4 * ,*4  4 •«  f-  Ml  - 41X4.4*4*44*  — 4 — — 4 — — 4 0 — 40—40—  4 

oooeoooooooocoooooocoooooooooocoocococeocooooooocco 
• ••  •••«•••••  •• 
S'  o • a •*  •>*#  c mmmujmm  mmm  wmmmmm  u*u*ui  £*£***  uim  m m m m ww  uj  wm  uj  In  m In  * it 


► u 

— u 

OU  «\ 

MW  */»  W> 

- V ft 

* — a*i  4 


<rt  W 3 s u 
4 u m u • u 

— V * V *v 
* — UJ  * - 

fi 


o 

♦ $ 

* O 

4 C 

4 — 

4 — 

m — 

m -4 

m 44 

m — 4 

‘ N 

m mi 

4 Ml 

4 ix 

4 Ml 

4 Ml 

4 IX 

4 

o 

- o 

- O 

— o 

-4  O 

— O 

- o 

— o 

4-  O 

- O - 

> o — © 

— o 

— o 

- O 

— o 

— O 

1 

• 

• 

1 • 

• 1 

• i 

• • 

• • 

• 1 1 

• 

i • 

• • 

• • 

1 « 

• • 

1 1 

• 

Ul 

Ul  Ul 

Ul  u. 

Ul  Ul 

Ul  Ul 

Ui  Ui 

Ui  Ui 

Ui  Ul 

Ul  Ul  u 

Ul 

a a 

Ul  Ul 

Ul  Ul 

Ui  Ua 

Ul  u. 

Ui 

mi  — 

m 4 

O m 

4 IM 

m 4 m«  in  o • 

Ml  -* 

— M>  - 

4 9 

Ml  4 4 4 

m mi 

4 O 

4 

• 

• 

Ml 

o © 4 4 

• • • • 

4 IX  m — 

4 O m mo  mi 
• • • • • • 
mi  — 4 — m 

— M 
• • 
M IM 

4 4 

• • 

4 Mi 

4 4 

• • 

IM  44 

4 M M 

• • 1 

Ml  -4  - 

n» 

• 

4 

Ml  4 
• • 

— 4 

• 4 

• • 
4 4 

4 m 
• • 

4 m 

4 m 
• • 
4 Ml 

4 4 

• • 
IM  — 

m 4 
• • 
♦x  — 

» 

• 

4 

eS?ooo3So 

4 — m 

4 — 4 

4-4 

40  4 

40  4 

ooo 

4 O 

o 

oo  6 oo  6 ooooob 

oo  b bo  b o 

oo 

oo  o o o o 

ooo  oooooo 

o O 

• 

• 

• 

• 

• 

1 

1 

• 

• i 

1 

• 

• 

» 

• 

• 

• 

hi 

Ul  Ul  Ul 

Ul  Ul  Ul 

Ul  Ul  Ul 

*** 

Ul  1*1 

4 0 4 

Ml  © Ml 

9 * * 

4 4- 4 

m+  - 

4 4 IM 

MOO ► 440 

• 4 

— 44 

4 40 

m 40  4 4 

4 

4 4 + 

***** 

OMO 

m 0 4 

o o n-  mo  *4 

4 4 Ml  M»  4 K 

mi  m o mi 

44 

4 4 M|  4m  4 

— m 

- 

f«  — r 

— — 4 

-44 

4 4 mi  m m *4 

IM  44  44 

mn>' 

— Ml  4 

— — m 4 4 — 

4 4 — 

4 Ml  4 

— M|  — 

O 

mo 

4 o 

* o 

m o 

m o 

m o 

mo 

mo 

mo 

m © 

m o 

m c 

m o 

m © 

m 

O 

o o 

o o 

o o 

o o 

o © 

o o 

o e 

oc 

o o 

o o 

o o 

o O 

o o 

o c 

o o 

o © 

o 

Ul 

o 

ss* 

4 o 

as 

Ui  t*l 

40 

?* 

ss 

tt 

UJ  Ul 
NO 

Ul  Ul 

oc 

Ul  41 

o o 

H 
• • 
•X  — 

Ul  u* 

4 O 

ui  ui 
o o 

Ui  Ui 
Ml  O 

Ul  Ul 

40 

Ui 

4 

o 

• 

♦ e 

• • 
• *4 

• o 

• • 
m — 

oo 
• • 

a© 

• • 

m o 
• • 
Mi  «• 

-4  O 
• • 
IM  44 

4 0 

• • 
m — 

P-  o 
• • 
IM  «4 

4 © 

• • 

IM  *4 

-o 
• • 
IM  -4 

• O 

• • 

IX  — 

— o 
• • 
Ml  — 

• o 

• • 

r o 
• • 

Ml  O 

• • 

O 

• 

IM 

4 m 

4 Ml 

m mi 

m -4 

m — 

m o 

mo 

me 

m — 

m — 

mix 

m mi 

m Mi 

m *m 

4 

c 

O C 

O O 

c o 

© o 

c c 

c c 

o o 

o o 

c c 

oc 

o c 

o o 

o e 

c c 

c c 

o o 

o 

• 

1 

• 

• 

• 

• 

I 

Ul 

Ul  Ul 

Ul  Ul 

Ul  Ul 

Ul  Ul 

Ul  Ul 

Ul  Ul 

Ul  Ul 

Ul  UJ 

Ul  Ul 

UlUa 

Ul  U! 

Ul  Ul 

Ul  Ul 

UJ  Ul 

Ua  Ui 

Ui  Ui 

Ui 

IT  a- 

4 « 

4.  4 

4 M 

m — 

•»  - 

4-  C 

Ml  C 

44  V 

c — 

!? 

4 m 

C 4 

IX  4 

c m 

4 

*x  4 

— m 

— M* 

40 

m 4 

44  4« 

4 Ml 

► r 

4 C 

— m 

4 4 

— 4 

4 4 

m #m 

44 

e 

IM 

M- 

m m 

— ix 

- 4 

v m 

IM  — 

m 4 

IM  Ml 

*M  IM 

m 4 

Ml  — 

X4 

Mi  m 

- - 

- 4 

— IX 

- 

4. 

— M 

44  IM 

o m 

O 4 

P 4 

— 4 

— 4 

— 4 

— 4 

Ml  4 

IX 

o 

o c 

o o 

o © 

O O 

o o 

o o 

O C 

o o 

o © 
l 

O O 

C C 

O O 

O © 

o o 

o e 

O O 

c 

1 

1 

T 

« 

I 

1 

• 

1 

1 

1 1 

• • 

• 1 

l • 

• • 

• 

UJ 

u u 

UJ  Ul 

Ul  u. 

u.  Ul 

Ul  Ul 

u.  u. 

U Ul 

Ul  Ul 

4 Ul 

UJ  Ui 

ui  Ui 

Ui  UJ 

U-  Ui 

UJ 

4 

O 4 

► o 

it 

4 4 

IM  M 

m m 

m m 

40 

r 4 

44 

4 4 

4m 

4 4 

« — 

— c 

44  4 

• 

ix 

4 a 

M ► 

4 Ml 

C 44 

• K 

4 |M 

Ml  4 

*m  m 

-4 

4 4 

4m 

m m. 

Mi  4 

Ml  4 

4 4 

4 

V 

«M4 

4 Ml 

4 r 

• m 

4- 

IM  |M 

Ml  — 

- 4 

4- 

4* 

4 4 

4 Ml 

IX  4 

— IX 

4 4 

4 

4 

44 

m 4 

m 4 

•M  4 

m 4 

m 4 

4 4 

im  y> 

*m  r 

4 4 

4 «r 

IM  4 

m 4 

m 4 

IM  4 

4 4 

4 

o 

c c 

C C 

o o 

O O 

c o 

o o 

O O 

O O 

o c 

c e 

o c 

o o 

C 6 

c c 

o o 

O O 

C 

u 

u.  u 

*$ 

Ul  Ui 

'tt 

Ui  u. 

Ul  Ua 

Ui  w 

Ul  Ui 

U.  Ua 

Ua  U. 

U.  Ul 

Ml 

mi  m 

m ir 

4 « 

4 ni 

4 m 

4 M 

4 4 

4 4 

4 4 

4 — 

4 4 

4 

« 

o c 

9 9 

O 4 

o c 

O 4 

O Ml 

O 4 

O — 

G im 

o - 

O 4 

©*X 

o o 

O 4 

e m 

C 4 

3 

4.  «4 

— 4 

— 4 

44  Ml 

- - 

^ "* 

44  M. 

— m 

44  IM 

— Ml 

~~ 

- 4 

— 4 

-*  m 

— '4 

- 

© 

— ■ 

s c 

— « 

c e 

Ml  • 

o o 

Ml  4 

o o 

M|  • 

o c 

Ml  4 

o o 

Ml  C 

o o 

Ml  M 

© C 

M,  M 

O O 

Ml 

c o 

1X4 
O O 

Ml  4 
C C 

IX  4 

o e 

Ml  4 
C O 

Ml  4 

o e 

o o 

Ml 

O 

Ul 

u.  u< 

Ul  Ua 

Ul  UJ 

it 

Ul  1 *1 

Ul  ul 

Ua  ly 

Ul  Ul 

Ul  Ul 

Ui  Ui 

U.  Ui 

Ui  Ui 

Ul 

Mi 

•V  9 

4 — 

c - 

4 * 

• 4 

m 4 

« c 

* in 

4 © 

O 4 

m 4 

4 m 

« — 

4 — 

4 

— 

► « 

— 4 

*>  c 

4 * 

► M> 

4 r 

*M  C 

m 9 

t*  V 

m m. 

— C 

4 Ml 

— 4 

4 »x 

4 4 

9 m 

4 

- 

4 4 

a m 

- 4 

Ml  IX 

4 •* 

m -4 

m ► 

4 m 

4 4 

4 4 

IM  Ml 

m — 

Ml  mm 

IX  4 

— 4 

- 

4 

9 * 

9 *> 

C 4 

Ml  4 

♦x  4 

ix  4 

*x  m 

rx  m 

Ml  4 

o 

o © 

o o 

o © 

OO 

o o 

o o 

o o 

o o 

o o 

O O 

O O 

o c 

o O 

o o 

o o 

o o 

o 

1 

• 

1 

l 

• 

• 

• 

• 

i. 

1 

1 

1 

1 

I 

Ul 

Ul  Ul 

uj* 

Ul  Ul 

Ul  Ul 

Ul  Ui 

Ul  ul 

UIUj 

Ul  Ui 

Ul  Ul 

41  U. 

Ui  Ui 

uut 

Ul  Ul 

Ul  Ul 

ass 

Ui  Ul 

Ul 

4 

»n  * 

m * 

• — 

4 4 

« m 

4 4 

4 Mi 

« 4 

4 4- 

C 4 

4m 

C G 

m 4 

4 4 

o 

* 

— 4 

O M 

m «n 

-4  IM 

mi  n» 

m -4 

o • 

— m 

>4  **< 

4 4 

4© 

4-  44 

O 4 

IX  IX 

4 4 

Ml 

♦ 

m rx 

M M 

•4  4 

4 4 

4 44 

4 4 

m 4 

M,  IX 

SZ 

— 4 

4 4 

m mi 

4 — 

m m 

ix  n 

*4  — 

— 

• 

— • 

C • 

O 4 

O * 

C M- 

o r» 

4*  n* 

Mim 

mi  m 

mi  m 

ix  m 

Ml  4 

Ml 

o 

o o 

o c 

o o 

o o 

O © 

© o 

o o 

o o 

o o 

O O 

O O 

c o 

o o 

o o 

o e 

o o 

o 

• 

• 

• 

1 

1 

1 

1 

a 

• 

• 

i 

u. 

Ul  Ul 

ui  Ui 

Ul  u. 

«*•  ul 

Ul  UJ 

ui  ui 

Ui  ul 

Ua  U> 

M.  Ua 

Ul  UJ 

Ul  ul 

Ul  u. 

Ui  i*i 

St* 

Ua  Ul 

hf 

4 

4 * 

♦ o 

Ml  -4 

m in 

4 4 

Ml  —4 

O JM 

— m 

n m- 

4 4 

O M 1 

O 4 

m o 

M 4 

m 

r 

• 

4 Ml 

• • 

* • 

• • 

4 O 
• • 

-•  4 
• • 

m -4 

• • 

44  « 

• • 

• c 

• • 

Ml  • 

4 m 

44 

Ml  4* 

• • 

4 4 
• • 

4 4 
• * 

n ix 

• • 

m mi 

• • 

4m 

in 

— im 

e — 

4 44 

m m 

p4  m 

44  44 

4 44 

4 m 

M Ml 

4-  -1 

4*  44 

4 4 

4 m 

IM  Mi 

M|  — 

— 4 

— 

o 

o 

c 

©• 

c 

9 

o 

O 

o 

C 

O 

C 

c 

O 

O 

© 

o 

o 

o 

o 

o 

O 

o 

O 

o 

e 

o 

o 

o 

o 

o 

O 

o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

• 

• 

• 

• 

• 

• 

o 

o 

o 

© 

o 

O 

o 

c 

o 

o 

o 

o 

o 

o 

o 

o 

o 

4 

• 

c 

Mi 

i 

9 

e 

4 

4 

o 

IX 

4 

4 

• 

Ml 

M 

4 

#M 

IM 

m 

4 

4 

4 

4 

4 

m 

m 

m 

m 

m 

Table  1-3  (Continued) 


i 


\ 


S 


l 

J 


ft 

u 

5o^o 

ooooooc 
• I 

55 

000060060060066 
< 1 1 1 1 

6 — 06  — 06  — 00  —6  6 — 66m66  — 66  — 
1 1 1 1 1 • ii 

KJ 

ft  . 

2* 

Ui  uiuiui 

UJ  UJ  UI  UIUJ 
ft  ft  ft  ft  ft 

Ui  V 

I o 

o ■#  *r  ft 

ft  ft  ft 

ft  ft  O ft 

O M 

ft  M ft  ft  M ft 

ft  M ft 

ft  ft  ft 

4 ft  — 4 

M M ft  ft  ft 

ft  ftft  ft  ~4  ft  ► ft 
OOOO-ftOMO 

I V 

rs 

UJ  X 

e — 

ft  e ft  * 

Oft  M 
• • • 

OMM  O 

• • • « 

' ft  — 
> • • 

O —O  O ft  ft 

O ft  ft 
• • • 

• • • 

Oft  ft 

til 

O ft  M O 
• • • • 

O ft  OO  M 

0 — * Oft 

— - ft  — 

— ft  — 

M ft  M m* 

i mm  mrn 

mm  rnm  rnm  rnm  9 9 

— ft  ft 

— N ft  — M 4 

-•ft  - 

►- 

— — x — 

— M — 

— M x 

>4  M 

M*  ftM* 

X 

o u 

ocooooooooooooooooooooooo 
• II  II  II  II  II  II  II  11 

oeooeooeooeoeoo 
11  11  ii  ii  ii 

00000000 
• 1 1 7 • 

u 

V 

1“ 

Ui  UJUJUi 

Ui  Ui  Ui 
ft  ft  0 

ft 

ft  — ft  o 

♦ ir  0 - 

ft  O ftft  ft 
— • O — ft 

s 

ft  X 

ft  ft-  - 

O-ft 

M — M ft 

ft  -ft 

O -*  «•  ft  — f-  O 

M M • M M ft 

ft  X 

ft  Mft  ft  rnm 

ft  ft 

_JU» 

• • • • 

• • • 

• • • • 

• • • 

• • • 

• • • 

» t t 

• • • • 

• t • 

s 

ft  u 

• ft-  ft  *•  - — 

► — ♦ ft 

M«X— M*>  — — 

ft  — ft 

m mm  9 mm  mm 

» — ft  - 

M M — — M 

mm  mm  9 mm  mm  9 

MMM 

mm  M 9 — — 

► u 

ft  ft  ft  ft  ft 

X 94 

ft 

ft  ft  ft  ft 

ft  ft 

► ft 

► ft 

► ft 

4 • 4 

m ft  • ft 

ft  ft 

ft  ftft 

9 - o 

- O 

-«  O - 

O 

- 0 - 0 - 0 

mrn  0 

— 0 - 

*-  Vi 

1 1 1 

1 I 

• II 

• 

1 II  I 

1 I 

1 1 

1 • 

• 1 • 

• I I 

• li  • 

1 1 

1 1 I 

o u 

mm  X 

Ui  Ui  i* 

Ui  Ui 

U Ui  Ui 

Ui 

Ui  U Ui  UJ 

Ui  Ui 

Ui  Ui 

UJ  UJ 

UJ  Ui  Ui 

Ui  Ui  Ui 

UJ  Ui  Ui  Ui 

m w 

- - ft  ft  ft 

O Oft 

ft 

X ft  — ft  M ft 

M ft  ft  ft 

ft  ft  ft 

«*  "V 

ft  X 

ft  «♦  ft  o X 

O ft  o 

ft  ft  1-  ft  M 

• ft 

ftft  ft 

O M « ft 

• 0 

O ftft 
• • • 
ft  ftft 

r — 

*5 

• • • 

ft  - ft 

• • 
— ft 

• « < 
• ft  ft 

• 

M 

• • • 9 

9 mm  m mm 

• • 

M ft 

• • 
M ft 

• • 
— 4 

• • • 
— ft  — 

• • • 

Mft  M 

• • • • 

• —ftft 

• • 
ft  ft 

ut 

ooojoooobooooobooooooooQOooo 

ft  ft  M ft 
OOOO 

m mm  ^ m mm 

OO  OOO 

X ft  M x ft  M 

OOOOOO 

?$s 

SS8SS 

ft  X 

1 I 

1 1 

• 1 • 

1 

• I II 

• I 

1 • 

• I 

7 • • 

7 ii 

7 • 11 

II  II 

* V 

u 

Ui  Ui  UI  Ui 

Ui  Ui  Ui 

Ui  Ui  Ui  Ui 

U*  UJ  Ui 

Ui  UI  UJ 

UJ  U Ui  UJ 

O 3 ft  M ft 

Ui  Ui  Ui 

M m X-  — 9 

ft  ft 

ft  ft  ft  • ft  m 

ft  — O 

ft  ft  M 

ft  mm  m mrn  X M mm 

«•  X 

ft  X 

*■  M«  ft 

ft  O — O ft  ft 

ft  ft  ft 

— M - 

z — 

UJ  ft 

is 

”* 

• • • i 

ft  ft  ft  ft 

• • • 
49  4 

• • t • 
ft  ft  M ft 

• • 

rnm  mm 

ft  mm  ft  ft  ft 

• • • 
Mft  ft 

• II 
MM  ft 

9 9 9 
mm  m*  M 

• • • • 
— • — — 

M M ft  — — M 

• • • 
9 4 mm 

O 4 0 

ft  c 

ft  O 

ft  O 

ft  C ft  0 

S g 

ft  O 

ft  O 

ft  0 

ftp  ft  0 

ft  O ft  O 

ftp 

ft  O ft 

o o c 

o c 

o e 

o o 

CO  00 

0 e 

O O 

0 0 

00  0 c 

OO  OO 

OO 

e 0 0 

o sj 

— UI 

— X 

u 
♦ o 

O X 

Ui  Ui  Ui 
O MO 
O O O 

ui  ui 

sg 

Ui  Ui 

Sg 

UJ  UJ 

ft  o 
ft  o 

UJ  UJ  UJ  UJ 

ft  0 ft  0 

*S 

jg 

SS 

Ui  u» 

ft  0 

• 0 

UJ  Ui  UJ  Ui 

*g  28 

82  SS 

— O MO 

«S  SS  s 

x 5 ftp  ft 

• • • 

• • 

• • 

• • 

• 9 

• • 

— ft  — 

X mm 

ft  — 

ft  — 

ft  - ft  — 

ft  mm 

M mm 

M — 

rnm  mrn 

MM  MM 

mm  9 mrn 

* — 

• — • 

ft  4 4 

o o c 

ft  ft 

oc 

ft  ft 

o c 

SS 

SS  SS 

gg 

gg 

SS 

ft  ft 

0 0 

SS  SS 

SS  08 

SS 

ftp  ft 

00  0 

1 1 

1 

• 

1 

• 1 

1 

• 

1 

1 

1 1 

• 1 

1 

1 

M V 

Ui  UI  UI 

2* 

Ui  UI 

Ui  Ui 

Ui  Ui  UJ  Ui 

UJ  UJ 

UI  Ui 

Ui  Ui 

Ui  UJ 

Ui  Ui  Ui  UJ 

UI  Ui  Ui  Ui 

Ui  UI 

SS  s 

c o 

o 

O M ft 

ft  — 

iz 

« © ft  M 

ft  — 

ft  mm 

M 9 

ft  ft 

ft  — — - 

C — ft*- 

ft  ft 

U -V 

UJ  X 

M CM 

ft  * 

ft  c 

ft  O ft  9 

— 0 

ft  « 

M •- 

ft  ft 

ft  M —ft 

- ft  Mft 

X X- 

ftft  ft 

• • • 

• • 

9 • 

• • 

• • 

9 9 

• • • 

- ft  ft 

M M 

ft  - 

ft  ft 

ft  M ft  • 

ft  ft 

M — 

M ft 

Mft 

MM  — ► 

Mft  9 mm 

K K 

ftft  ft 

o 

ft  M ft 

M O 

ft  o 

ft  — 

ft  — ft  M 

ft  M 

ft  ft 

4 ft 

ft  ft 

ftft  ftft 

ftft  ftft 

4 ft 

ftft  ft 

5> 

O CO 

o - 

o — 

e — 

O - O - 

O - 

O - 

0 — 0 - 

1 1 1 

1 I 

• 1 

• • 

• 1 II 

1 1 

1 1 

1 1 

1 1 

II  II 

• 1 II 

1 I 

11  7 

m u 

u.  u.  ui 

U Ui 

U Ui 

U Ui 

WO  Ui  U 

SS 

UI  Ui 

Ui  U 

Ui  UI 

Ui  U.  Ui  u 

U.'  Ui  UJ  Ui 

Ui  UJ 

UJ  UJ  UJ 
ftft  ft 

o t 

V-/ 

M X ~ 

ft  ft 

ft  ft 

— ft 

mm  9 -■ 

• ft 

ft  ft 

ft  -■ 

mm  0 MM 

MM  • - 

ft  ft 

* X 

o o 

M ft  — 

ft  m 

— ft 

e ft 

C ft  ft  ft 

• ft 

M ft 

ft  M 

— ft 

ft  b-  C ft 

- M ft  - 

0 m 

ftft  — 

X 

M 

ft  X — 

— 4* 

ft  — 

« ft 

ft  - M ft 

mrn  mm 

— ft 

ft  M 

ft  ft 

ft  M ftft 

1 

1 

€ 

Z 

M ft 

4 mm  9 

ft  ft  ft 

o o o 

ft  ft 

c e 

ft  ft 

o o 

ft  ft 

o c* 

ft  ft  ft  Mi 

00  00 

ft  M 
O C 

ft  M 
O O 

ft  N 
O O 

ft  M 
OO 

ft  — ft  M 

00  0 e 

ft  mm  ft  — 

0 e 00 

ft  M 

O O 

ft  M ft 

c c e 

_ V/ 

ft  * 

u.  u.  u 

UJ  UI 

Ui  Ui 

SS  s 
• • • 

C U 

Ui  O 

x « C 

ft  ft 

4 ft 

ft  M 

ft  ft  ft  M 

ft  • 

ft  ft 

4 ft 

ft  ft 

Si  gi 

99  99 

4 4 49 

ft  9 

O Oft 
• • • 

O ft 

• • 

O ft 
• • 

o ft 
• • 

0 M O mm 
9 9 9 9 

O M 

• • 

O ft 
• • 

O M 
• • 

e ft 
• • 

0 9 Oft 
• • 99 

C ft 
• • 

Mi 

— - ft 

— ft 

— M 

mm  mm  rnm  m 

— ft 

rnm  m 

-•  M 

mm  M 

-ft  - ft 

Mft  m n# 

M M 

MM  M 

ft  M ft 

M ft 

— ft 

— ft  — ft 

— ft 

mm  ft 

Mft  Mft 

o o c 

o c 

e c 

C 9 

00  c c 

O O 

0 e 

O O 

O O 

00  0 e 

CO  OC 

p 0 

OO  O 

u 

o 

— Ui 

w ft1 

Ui  U.  Ui 
ft  ft  ft 

Ui  Ui 

ft  ft 

UI  Ui 

— 9 

UJ  U 
• M 

Ui  Ui  Ui  Ui 

IS  St 

Ui  tu 

« - 

Ui  Ui 

C ft 

UJ  U 
ft  ft 

Ui  Ui 

r • 

Ui  Ui  Ui  Ui 

ftft  ftft 

Ui  Ui  Ui  Ui 
0 C Mft 

UJ  Ui 

■ M 

SS  s 

O X 

— * 

• ft 

M ft 

9-  mm 

M ft 

ft  J 

ft  ft 

ftft  — ft 

C — ftft 

OX  M 

• • • 

9 9 

• • • • 

9 9 

• • 

ft  — ft 

— M 

ft  - 

9 mm 

ft  - ft  ft 

ft  ft 

ft  ft 

ft  ft 

MM 

M — Mm 

mm  mrn  M « 

- ft 

— ft  ft 

M ft  M 

ft  M 

ft  M 

ft  ■-> 

ft  mm  ft  mm 

ft  — 

ft  O 

ft  O 

O O o 

O C 

O O 

o O 

OO  OO 

O O 

O O 

O O 

O O 

OO  OO 

OO  OO 

P O 

OQ  O 

_ W 

1 

• 

1 

• 

• 1 

1 

1 

1 

1 

1 It 

• 1 II 

1 1 

• 7 • 

V 

r uj 

UI  Ui  Ui 

U U 

Ui  Ui 

Ui  Ui  UI  UI 

uj  u« 

UJ  UJ 

UJ  Ui 

Ui  Ui 

UIUI  Ui  y 

SS 

5f w 

o x 

-»  X 

« ft  ft 
ft  ft  M 

ft  ft 
ft  c 

4 ft 
« ft 

ft  ft 
ft  f* 

ft  « ft  M 
ft  • ft  ft 

ft  M 

ft  ft 

M M 
— ft 

ft  ft 

M ft 

MC 

ftft  Oft 

• ft  Oft 

ft  M ft 

Mft  ft 

ftft  ftft 

ft  ft 

• • • 

• • 

• • 

• • 

• • • 

ft  ft  ft 

ft  M 

ft  — 

ft  ft 

M ft  — M 

Ml 

mm  ft 

ft  4 

ft  M 

ft  — ftp 

Mft  MM 

— ft  ft 

ft  ft  ft 

ft  ft 

ft  ft 

ft  ft 

ft  ft  ft  ft 

ft  ft 

ft  *• 

4 M 

ft  — 

O CO 

o o 

o o 

o o 

00  O C 

O O 

O O 

O O 

O O 

00  00 

OO  OO 

p 0 

00  e 

• 

• 

• 

1 

• I 

1 

1 

1 

• 

• I 

1 I 

7 

• 1 

u 

5 Si 

U.  Ui  Ui 

u.  u 

w*  u 

Ui  UI  Ui  Ui 

UJ  UJ 

Urn  Ui 

Ui  UJ 

Ui  Ui  Ui  Ui 

— • ftft 

Ui  Ui  Ui  Ui 
ftft  M — 

Ui  Ui 

S t t 

ft  ft  M 

• ft 

ft  — 

- ft 

• • • ft 

ft  ft 

ft  — 

ft  ft 

e 0 

MM 

ft  — ft 

« ft 

M M 

— ft 

MM  ft  ft 

M — 

mrn  mm 

• M 

- « 

• ft  ftft 

- • ft  ft 

MM 

• • • 

ft  ft  M 

• • 
ft  — 

• • 
ft  ft 

• 9 
ft  ft 

• • • 1 

M ft  — - 

• • 
mm  mrn 

9 9 

9 ft 

• • 
ft  4 

• • 
ft  M 

t • • • 

ft  — Mft 

• • • • 
Mft  — ft 

• • 
m M 

• 9 9 

ft  — ft 

9 c* 

o 

C 

O C 

c 

c 

C 

C 

O P 

9 O 

O 

© o 

o 

o 

© 

O O 

0 

O 

O 

8 

© 6 

0 5 

P 

O 

• • 

• 

• 

• 

• • 

• 

• 

• 

9 o 

O 

o 

o 

O O 

0 

O 

C 

P 

0 0 

0 0 

P 

O 

© M 

ft 

ft 

ft 

O M 

ft 

ft 

ft 

O 

M 4 

• • 

O 

M 

ft  ft 

ft 

ft 

ft 

ft  ft 

ft 

ft 

• 

• • 

• • 

ft 

ft 

23 


Table  1-3  (Continued) 


l 

t 


¥8 


Z 

8 


n 


O O—OO  — O O - C O — 0 0-00-  oo-oo—  o o — 


a 4 w wi 

So  x a x ■. 

a- 


»«o«k«  4A««4X4—  x 4 4X 


*i 


OOOOOOOOOOOOO 

• 1 1 1 1 1 • 1 i 

4 — — X — — — — — A X — • • 


» 4 x x o * 
>0  00  — c 
l • • 

HI  HI  Hi  HI  HI  H 

« OH  a — 0 


A 

• A 4 #| 

• A> 

• 4 • 4 

» • 40  X 4 X 

4 

► u 

0 

- O 

- O 

— O 

• — O 

- O 

— O — 

*-  0 

• 

• • 

1 • 

• • 

• • • 1 

t I I 

• I 

• • 

• 

0 w 

— N. 

Hi 

Hi  Hi  Hi  Hi 

Hi  Hi 

HJ  HJ  Hi  Hi 

S 88  8 

8 8 

Hi 

NU 

A A 

O 

O — 

• • 

IS*  AS 

X 4 4 4 
• A A f 

X 

— ■ x 

Z % 

O 

4 IS* 

4 A 

wm  wm 

4 As 

4 4 

4 4 

X 

Z - 

*5 

• 

A 

• • 
a x 

• • 
A *4 

0 0 

AS  p« 

• • • 

X AS  X X 

» • • 

• • 
— X 

• • 
-A  A 

0 

0 0 4 0 

A40AA0AA—  A A- 

000000000  000 

> • As  — 

MX 

4 4 X 

4 X 

HI 

O 

OOO 

OOO 

OOO 

OOOOO 

« X 

• 1 

• I 

1 1 

• • • • • 1 

• 1 I 

1 1 I 

• 1 I 

• • 

asj 

31 

U 

U! 

Ski  Ski  HJ 

Hi  Hi  Hi 

Ski  Hi 

♦ 0 

A u 

0 u 

4 AS  m IS*  — 4 x x 4 

1 SIS  0 AS  A h-  0 

4 4- 

— 4 

— V. 

A X 

t x 

4 

4 — X 

A Ah 

4 4 A 

* - 

HiZ 

• 

• • • 
0*S|  # 

• • • 

in  r*  # 

1 • • • 

A ru  jj? 

0 • • 

_•  • 

t O 

- 

< 4- 

1 

I • 

O 

A O 

AS  0 

ASO 

AS  © AS  O AS  O 

xo 

XO 

X 

O 

O O 

0 c 

O O 

00  00  0 c 

0 0 

00 

O 

a <j 

w 

HJ 

Hi  HI 

¥8 

Hi  Hi 

Hi  HI  HI  Hi 

1 HI  Hi 

Hi  Hi 

Hi  Hi 

Hi 

O 

• O 

ISI  O 

00  CO 

4 O 

A O 

• 

— X 

b x 

O 

US  O 

• e 

X 0 

X O 4 0 

A O 

40 

4 

0 - 

• 

• • 

• • 

• • 

• • • l 

t • • 

• 0 

• 

♦ - 

AS  -4 

•A  ,4 

X mm  mm 

AA  — 

4 — 

4 — 

4 

4 

♦"  C 

a 0 

ASO 

AS  rnm  AS  x 

1 A X 

X A 

X A 

X 

C 

c — 

0 - 

O - 

O — O - 

, O - 

O — 

O — 

O 

• 

1 

• 

• 

• ( 

• 

1 

X U 

V 

Hi 

Hi  HJ 

Hi  Hi 

Hi  HJ 

W Hi  WH 

H Hi 

Hi  Hi 

Hi  HJ 

Hi 

O \J 

0 

« IN* 

4 4 

X C 

• 4 4 

4 4 

AO 

« 

u >. 

H*  >. 

X 

IT  O 

4 X 

X 4 

X — 4 4 

A X 

4 A 

4 A 

4 

m* 

• 

0 • 

• • 

• • 

4 • 

a a 

AS  -■ 

X 4 — 4 

— X 

4AS 

4 — 

4 

z 

X 

A X 

AS  • 

As  A 

AS  4 40 

4 - 

4 X 

X A 

X 

Hi 

O — 

0 — 

O - 

0 — OX  ox 

C X 

O X 

C 

A 

1 

» • 

I 1 

• 1 

II  II 

1 1 

1 1 

1 I 

• 

XU 

8 

Hi  Hi 

H.  Hi 

Hi  H> 

iy  H<  ly 

1 Hi  Hi 

Hi  Hi 

Hi  Hi 

HI 

00 

U 

a a 

a x 

0 4 

— X 

4 4 

4 4 

4 

r v 

ou 

1 r 

■0  a 

IT  4 

• « 

XX  X 4 

— X 

A 4 

4 A 

X 

• 

• 0 

• • 

• • 

• • 

0 0 

• 

X 

— 

4 

a mm 

X 4 

— • -4 

— X 4 X 

A A 

— A 

X 4 

A 

O 

a 0 

AS  C 

ASC 

A«  — AS  mm 

A mm 

AX 

AX 

A 

O 

0 c 

O O 

e 0 

00  c c 

c 0 

c c 

CO 

O 

ft  4 

1 1 

1 

• 

1 

u 

1 

Hi 

H.  HJ 

Hi  *j 

U.  H* 

u.  H.  H>  Ha 
4 4 4 - 

H.  Ski 

Iki 

5 u 

HiO 

A 

4 4 

4 X 

4 4 

4 4 

f X 

4 A 

4 

2 x 

*-  Hi 

C 4 

0 C 

O 4 

C X O As 

0 X 

O 4 

0 

a 

• 

• • 

• • 

• • 

* * • 1 

• • 

• • 

• • 

0 

Hi 

** 

— 4 

p-  AS 

"*  — 

— m — as 

1 - 

- 4 

— X 

♦ 

O 4 

O 4 

C 4 

O AS  O A 

C A 

O A 

O A 

O 

O O 

c 0 

O C 

CO  C C 

c 0 

O O 

O O 

0 

U 

1 

u 

— HJ 

Hi 

Hi  U, 

Hi  Hi 

Hi  Hi  Hi  Hi 

AS  « 4 - 

u 

uu  *r 

A 

C 4 

c c 

• — 

4 4 

X 4 

4 4 

A 

C v 

> x 

O 

► 1 

As  — 

4 4 

4 As  X A* 

A X 

4 4 

4 

• 

• • 

• • 

• • • 1 

• • 

• • 

• • 

• 

A 

x IM 

4 X 

AS  — 

4 4 A A* 

X A 

■*  — 

PA  — 

4 

IS* 

A M 

As  is* 

AS  X 

4 A AS  A 

4 4 

4 4 

4 A 

4 

O 

0 0 

O O 

O O 

O O O O 

OO 

O O 

O 

u 

1 

• 1 

• 1 

• • 

II  II 

• I 

1 1 

1 1 

• 

u 

* Hi 

Hi 

HI  HI 

Hi  HJ 

HI  HJ 

Hi  Hi  s«i  Hi 

Hi  Hi 

HI  Hi 

Hi  Hi 

Hi 

IS*  U 

— A 

a — 

c 0 

— AS 

— • 14 

4 X 

C A 

— 4 

O 

0 X 

> X 

A 

A — 

O • 

• O 

4 A 1-  X 

4 a 

4 A 

A 4 

• 

• • 

• • 

• 0 

• • • 1 

• • 

• • 

• • 

• 

«r 

4 A 

AS  — 

AS  — 

X A -a  mm 

4 4 

A — 

A A 

X 

0 

A0 

a© 

5 c 

AS  — A — 

4 x 

4 X 

4 

0 

O O 

0 5 

0 5 

OO  OO 

O O 

c 0 

O 

u 

1 

• 

• 

II  II 

• I 

1 I 

• 1 

• 

w 

IS*  W 

i* 

H* 

A 

Hi  H- 

O is* 

HI  HI 
AS  4 

(m  h • 
X x 

Hi  H*  Hi  H> 
4 X X O 

Hi  H. 

A 4 

H-  H> 
4 X 

Hi  ski 

Hi 

4 

Z X 

> X 

4 

A IS* 

O IS* 

0 0 

AO  A -a 

C X 

X 4 

HO 

4 

• 

• • 

• • 

• • 

• • • « 

• • 

0 0 

• 

O 

A A 

4 *A 

AX 

XX  rnm  AS 

4 -A 

4 4 

XX 

O 

O 

O 

O 

O 0 

O 

C 

C 

X 

O 

O 

O 

O 

O O 

O 

O 

0 

0 

• 

• 

• 

• « 

• 

4 

O 

O 

O 

P 

O O 

O 

O 

0 

f 

4 

4 

Q 

4 • 

4 

C 

4 

4 

4 

O 

— 

0 0 

•A 

X 

24 


REFERENCES 


1-1.  Hamlin,  D.A.  and  M.R  Schoonover,  The  ROSCOE  Manual , Volume  14a: 
Ambient  Atmosphere  (Major  and  Minor  Neutral  Species  and  Iono- 
sphere), DNA  3964F-14a,  Science  Applications,  Inc.,  13  June  1975. 

1-2.  Myers,  B.F.,  The  ROSCOE  Manual,  Volume  14b:  Midlatitude  Density 
Profile  of  Selected  Atmospheric  Species,  DNA  3964F-146,  Science 
Applications,  Inc.,  13  June  1975. 

1-3.  Hamlin,  D.A.,  et  al,  ROSCOE  Program  Listings  Obtained,  1 March 
1978. 

1-4.  Booker,  H.G.,  "Fitting  of  Multi-region  Ionospheric  Profiles  of 
Electron  Density  by  a Single  Analytic  Function  of  Height," 

«/'.  Atmos.  Terr.  Phys. , Vol.  39,  pp  619-623,  1977. 


SECTION  2 

AMBIENT  IONOSPHERE  AND  ION  NEUTRAL 
COLLISION  FREQUENCY  MODELS 


INTRODUCTION 

Simple  models  of  the  daytime  (noon)  and  nighttime  (midnight) 
normal  ionospheres  are  used  in  the  WEPH  code  to  indicate  when  nuclear 
burst  produced  ionization  has  fallen  to  normal  (ambient)  levels  and 
to  provide  a reference  for  electromagnetic  propagation  effects.  The 
models  have  generally  been  satisfactory  for  propagation  frequencies 
above  the  HF  band  where  propagation  effects  are  dependent  on  the  in- 
tegral of  ionization  along  the  propagation  path,  but  can  result  in 
significant  errors  at  lower  frequencies  where  the  gradient  of  ioniza- 
tion is  important  (effects  of  striations  and  local  inhomogeneities 
are  neglected  in  this  discussion). 

Effective,  D-region  ionosphere  models  have  recently  been  pre- 
pared that  provide  improved  predictions  for  propagation  in  the  VLF 
and  LF  bands  (References  2-1  and  2-2).  These  models  are  used  in  con- 
junction with  E-  and  F-region  ionosphere  models  developed  for  use  in 
HF  propagation  predictions  (Reference  2-3)  to  prepare  a new  ambient 
ionosphere  model  for  the  WEPH  code. 

A simple  model  for  the  ion-neutral  collision  frequency  has 
also  been  used  in  the  WEPH  code  due  to  uncertainties  in  theoretical 
and  experimental  data.  While  there  is  still  considerable  uncertainty 
in  specifying  the  ion-neutral  collision  frequency  the  model  has  been 
reformulated  in  terms  of  ion  mobility.  Mobility  values  are  used  which 
are  consistent  with  recent  measurements  and  which  result  in  ion- 
neutral  collision  frequencies  similar  to  those  used  at  the  Naval 
Ocean  Systems  Center  (NOSC)  to  obtain  agreement  with  ELF  propagation 
data. 
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CURRENT  WEPH  MODELS 

In  the  current  WEPH  code  the  ion-pair  production  rate  for 
altitudes  below  90  km  is  obtained  from  relations  developed  by  Science 
Application  Inc.  (SAI)  for  ROSCOE  (Reference  2-4).  The  electron  and 
ion  densities  for  ambient  conditions  are  found  from  the  ion-pair  pro- 
duction rate  and  the  I)- region  chemistry  model.  For  altitudes  above 
90  km  the  electron  density  is  obtained  from  relations  given  in  Refer- 
ence 2-4  and  an  effective  ion-pair  production  rate  derived  using  the 
E-  and  F-region  chemistry  model.  Ambient  models  for  he  electron 
and  nitrogen  vibrational  temperatures  arc  also  obtained  from  rela- 
tions given  in  Reference  2-4.  The  ambient  ionosphere  model  is  inde- 
pendent of  sunspot  number,  geographic  location,  season,  and  local 
time  (other  than  day  or  night). 

The  ion-neutral  collision  frequency  in  the  current  WEPH  code 
is  calculated  in  terms  of  the  average  electron-neutral  collision 
frequency 

v 

— em 

Vim  = 20  ’ 

For  altitudes  below  about  100  km  the  electron  neutral  collision  fre- 
quency can  be  expressed  in  terms  of  number  density  and  temperature 
or  pressure 

><>"• 

The  average  electron-neutral  collision  frequency  is  the  high  frequency 
average  (for  use  when  w >>  v)  defined  by  Shkorofsky  'Reference  2-5). 
The  average  ion-neutral  collision  frequency  is  the  low  frequency 
average  (for  use  when  co  <<  v) . 

D- REGION  IONOSPHERE  MODELS 
Theoretical  Models 

Theoretical  models  for  the  D-region  can  be  developed  by  formu- 
lating ion-pair  production  rates  and  atmospheric  chemistry  models. 
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In  general,  both  the  production  rate  and  chemistry  models  arc  depen- 
dent on  solar  conditions,  geographic  location,  season,  and  local  time. 
If  an  ambient  electron  density  profile  is  known  or  assumed,  it  can  be 
used  in  conjunction  with  an  atmospheric  chemistry  model  to  derive  an 
effective  ion-pair  production  rate.  This  will  insure  that  the  desired 
ambient  electron  density  is  obtained  when  the  ion-pair  production  rate 
is  equal  to  the  effective  value.  However,  if  the  chemistry  model  is 
in  error,  then  the  effective  ion-pair  production  rate  will  be  in 
error  and  calculations  of  the  transition  from  disturbed  to  ambient 
conditions  can  be  distorted. 

Figure  2-1  shows  an  example  of  daytime,  D-region  ionization 
rates.  Below  about  90  km  galatic  cosmic  rays  and  _yman  a radiation 
dominate.  Relations  for  these  sources  arc  given  in  the  literature 
(sec  for  example  References  2-6,  2-7,  and  2-8).  For  galactic  cosmic 
rays  the  daytime  and  nightime  ion-pair  production  rate  can  be  approxi- 
mated from  relations  given  in  Reference  2-8  as 


Figure  2-1.  Daytime  D-region  ionization  rates. 


28 


(2-1) 


qc  = q + R cos  x 


where 


x 


Y 

y 

R 

q 


'I  2n  Y < 4.5 

= < 

~ 2tt  4.5  < Y < 11 
= MOD  (y-9,11) 

= number  of  the  year  in  the  1900 
(eg,  1974  becomes  74) 

2^cmax  %min^ 

= i(q  + q . ) 

2 xmax  ncmm' 

/ 0. 056  + 0.7  sin40)  ab 


qcmax 


([(0.107  + 0.7  sin40)a  - o.8]b 
(0.  056  + 0,48  sin40)ab 


^cmin 


a 


([(0.12  + 0.48  sin40)a  - l]b 
4 x 10’^n 


a > 1 

a > 1 

a < 1 

a > 1 


b = exp(-6  x 10-6p) 

_2 

n = number  density  of  neutral  particles  (cm  ) 

. 2 

p = atmospheric  pressure  (dynes  cm  ) 

0 = minimum  (magnetic  latitude,  60°) 

Note  that  the  galactic  cosmic  source  is  maximum  at  sunspot  minimum 
and  minimum  at  sunspot  maximum. 


29 


For  Lyman  a radiation  the  daytime  ion-pair  production  rate  is 
given  in  Reference  2-7  as 

[,  . J] 

\a  =\ 

V^La^max  [*  * “74  ] 

where 

(qLct)max  = 6 x 10- 7 [NO]  exp(-  0.047p  secX) 

X * solar  zenith  angle 

_2 

p = atmospheric  pressure  (dynes  cm  ) 

[NO]  = nitric  oxide  number  density  (cm  ) 

At  night  scattered  Lyman  a radiation  may  be  important,  but  the  decay 
of  daytime  ionization  probably  dominates  above  about  70  km. 

A relatively  detailed  D-region  chemistry  model  is  used  in  the 
WEPH  code  to  determine  electron  and  ion  densities  for  a given  ion- 
pair  production  rate.  Electron  and  ion  densities  are  dependent  on 
the  presence  of  minor  neutral  species  (eg,  0,  02(A),  NO,  NO2,  I^O) 
as  well  as  the  major  species.  However,  the  model  was  developed  for 
disturbed  conditions  and  may  not  include  sufficient  detail  to  model 
the  atmospheric  response  to  low-level  ionization  sources.  Further, 
the  current  minor  neutral  species  model  (obtained  from  Reference  2-4) 
does  not  include  variations  with  solar  activity,  geographic  location, 
season,  or  local  time  other  than  day  or  night. 

Empirical  Models  for  Low 
Frequency  Propagation 

Davis  and  Berry  (Reference  2-1)  have  analyzed  over  600  elec- 
tron density  profiles  published  in  the  literature  and  have  derived 
an  electron  density  model  for  use  in  predicting  VLF  and  LF  propaga- 
tion. The  model  is  based  on  the  assumption  that  VLF  and  LF  propagation 


0 < Y S 4 


4 < Y < 11 


(2-2) 
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in  the  ambient  ionosphere  can  be  predicted  with  an  exponential  elec- 
tron density  profile  of  the  form 

Nc  = Nw  cx^aw(h  " hw}1  (2'3) 


where 

N = 1.43  x 107  exp (-  0.15h  ) 
w 1 w 


and  h and  a are  parameters  chosen  to  fit  the  electron  density  pro- 
w w 

files  in  the  altitude  region  where  VLF  and  LF  reflection  occurs.  This 
region  is  located  near  altitude  hQ  defined  by 


3,65 

w a +0.15 
w 


(2-4) 


The  expressions  for  h and  a given  in  Reference  2-1  are 

w w 

h = 71.81  - 7.84X.  + 8.04X.  - 1.23X,  - 
w 12  3 

0.0371X4  - 7.03XS  (2-5) 

a = 0.353  - 0.120X,  - 0.072X.  + 0.171X  (2-6) 

W i j J 

where 

Xj  = cos  X,  the  solar  zenith  angle 
X^  = cos  9,  the  geographic  latitude 

X^  = cos  <j>,  a seasonal  variable,  <p  = m ■ yy-1  3(2tt) 
m = month  number 

X^  = SSN,  the  Zurich  smoothed,  relative  sunspot  number 
for  month  of  prediction 

Xj.  = absorption  index,  with  values  0.0  for  quiet  condi- 
tions, 1.0  for  disturbed  conditions. 

Table  2-1  shows  predicted  values  for  h and  a given  in  Reference 

w w 

2-1  for  selected  conditions. 
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Table  2-1.  Predicted  values  of  and  from  Equations  2-1  and  2-2 
for  selected  conditions  (quiet  magnetic  absorption  index). 
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In  deriving  Equations  2-5  and  2-6  from  the  electron  density 
data  a propagation  frequency  of  50  kHz  was  used  to  define  the  reflec- 
tion altitude  hQ.  Use  of  a lower  or  higher  frequency  would  change  hQ 
and  thus  the  portion  of  the  measure  electron  density  profile  modeled 
by  Equation  2-3.  While  for  typical  electron  density  profiles  the  change 
in  hQ  would  be  only  a few  kilometers,  the  electron  density  gradient 
could  change  significantly.  Only  a few  of  the  profiles  used  to  de- 
rive Equations  2-5  and  2-6  were  for  nighttime  conditions.  Also,  at 
night  the  reflection  altitude  is  high  and  the  electron  density  re- 
quired to  cause  reflection  is  relatively  low.  Thus  measured  data 
near  the  reflection  altitude  can  have  considerable  uncertainty. 

Finally,  the  variation  of  electron  density  during  sunrise  and  sunset 
can  be  considerably  more  rapid  than  indicated  in  Equations  2-3,  2-5, 
and  2-6,  particularly  bolow  about  70  km  where  changes  in  minor  neu- 
tral species  can  affect  deionization  chemistry. 

Morfitt  (Reference  2-2)  has  derived  effective  exponential 
electron  density  profiles  from  analysis  of  VLF  and  LF  field  strength 
measurements.  The  field  strength  data  include  several  radials  out- 
ward from  Hawaii,  several  transmission  paths  across  the  continental 
United  States,  and  high  latitude  propagation  over  the  Greenland  Ice 
Cap.  Table  2-2  shows  effective  values  of  hw  and  given  in  Refer- 
ence 2-2  for  daytime  conditions.  Also  shown  are  the  values  of  h and 

w 

aw  predicted  by  Equations  2-5  and  2-6  for  similar  conditions.  Table 
2-3  shows  similar  data  for  nighttime  conditions  (nighttime  data  were 
only  taken  for  the  Pacific  paths).  The  value  of  has  a relatively 
strong  dependence  on  frequency  suggesting  that  the  actual  slope  of 
the  electron  density  profile  increases  with  altitude  in  the  reflec- 
tion region. 

E-  AND  F-REGION  IONOSPHERE  MODELS 
Theoretical  Models 

Theoretical  models  for  the  E-  and  F-regions  can  also  be  de- 
veloped (see,  for  example.  Reference  2-9),  but  the  formulation  is 


\ 
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Table  2-2.  Effective  values  of  hw  and  for  daytime  obtained 
from  VLF  and  LF  propagation  measurements. 


Propagation 

U 

Path  Location 

Season 

hw 

Midlatitude- 

Summer 

70 

Pacific 

Winter 

74 

Mldlatltude- 

Summer 

72 

U.S. 

Winter 

72 

Greenland 

Summer 

72 

Ice  Cap 


(Equations  2-3  and  2-4) 


“w 

hw 

“w 

0.35 

71.5 

0.3 

0.15 

72.4 

0.22 

0.15 

71.3 

0.31 

0.15 

72.9 

0.23 

0.15 

69.5 

0.33 

Table  2-3.  Effective  values  of  hw  and  for  nighttime  obtained 
from  VLF  and  LF  propagation  measurements.  Midlatitude 
(Pacific),  winter. 


Frequency 

(kHz) 

hw 

“w 

Below  10 

87 

0.15 

10-15 

87 

0.25 

15-25 

87 

0.35 

25-30 

88 

0.45 

30-40 

88 

0.55 

40-60 

88 

0.65 

(Equations  2-3  and  2-4) 


w 


a. 


w 


81.3  0.35 
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more  complex  and  is  dependent  on  phenomena  such  as  diffusion  which 
arc  difficult  to  model. 


Empirical  Models  for  HF  Propagation 

Stanford  Research  Institute  (SRI)  prepared  E-  and  F-region 
electron  density  models  for  use  in  the  NUCOM  code  that  are  based  on 
Institute  for  Telecommunication  Sciences  (ITS)  data  obtained  from 
vertical-incidence  sounders  (Reference  2-3). 


At  altitudes  between  90  km  and  the  peak  of  the  region  the 
electron  density  is  specified  by  three  profiles  of  the  form 


y > h > h - y (2-7) 


m 


where  h is  the  height  of  the  profile  maximum,  N is  the  electron 
density  at  that  maximum  and  ym  is  the  profile  semithickness.  Two  of 
the  profiles  model  the  E-region  and  the  F^  region  and  values  for 

Nmax’  hm’  antl  ym  are  obtainec*  ^roni  *TS  data  for  tbese  regions.  A 
third  profile  is  used  as  a filler  or  connecting  region  and  the  values 

for  N , h and  y are  derived  in  terms  of  the  parameters  for  the 
max  m ' m 1 

E-  and  F^-regions.  For  altitudes  above  the  peak  of  the  F^-region  the 
electron  density  is  computed  from 


N (h)  = N 
e 


max 


h - h / h - h 

<nwt ' o^nr)Jj  <2-8' 


h > h 

m 


where  N . h , and  y are  the  values  for  the  F. -region, 
max  mm  2 

For  the  E-region  the  parameters  are  taken  as 

h = 115  km 


max 


y 


m 


N 


max 


= 25  km 

= 1.24  x 104fJ 
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where  f_„  is  the  E- region  critical  frequency  obtained  from  the  ITS 
Ufc 

data.  For  the  F--region  the  parameters  h , y , and  N are  all 
obtained  from  ITS  data.  The  filler  or  connecting  region  (called  the 
Fj-region)  is  derived  by  requiring  that  it  overlays  with  the  F^-region 
by  half  the  semithickness  of  the  F^-region  and  by  choosing  the  bottom 
of  the  Fj -region  as  130  km.  The  peak  electron  density  of  the  Fj- 
region  is  related  to  the  peak  density  of  the  E-rcgion  by  requiring 
that  the  Fj -region  critical  frequency  be  twice  the  E-region  critical 
frequency  when  the  semithickness  of  the  F^-region  is  120  km.  The 
formula  used  is 


f ~ f QL 

r0Fl  0Ea 


(2-9) 


where 


n = 120 

and  y is  the  F -region  semithickness, 
m 1 & 

Tables  2-4,  2-5,  and  2-6  show  values  of  fQE>  f0p2»  hmF2’  and 
ymP2  obtained  from  ITS  data  for  selected  conditions. 


NEW  IONOSPHERE  MODEL 
General 

The  new  ionosphere  model  is  similar  to  the  current  WEPH  model 
in  that  effective  ion-pair  production  rates  are  derived  that  produce 
specified  electron  density  profiles  when  used  with  the  WEPH  code  chem- 
istry models.  For  the  E-  and  F-regions  electron  density  data  obtained 
from  the  SRI  model  are  used.  For  the  daytime  D-region  the  electron 
density  profile  obtained  from  Davis  and  Berry  (Equation  2-3)  is  used 
to  normalize  cosmic  ray  and  Lyman  a production  rates.  This  is  done 
since  only  a limited  altitude  region  is  described  by  Equation  2-3. 

For  the  nighttime  D-region  the  electron  density  obtained  from  Equa- 
tion 2-3  is  used  to  derive  an  exponential  production  rate  for  the 
upper  D-region  that  can  be  combined  with  the  cosmic  production  rate. 
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Table  2-4.  E-  and  F-region  data  for  March,  40  degrees  longitude. 


Sunspot  Number 

10 

3 

130 

Parameter 

Latitude 

Local  Time 

Local 

Time 

0 

12 

0 

12 

0 

0.288 

3.731 

0.451 

4.549 

f0E 

30 

0.177 

3.688 

0.318 

4.386 

(MHz) 

60 

0.393 

2.805 

0.557 

3.314 

0 

6.776 

8.729 

11.584 

12.756 

f0F2 

30 

3.611 

7.251 

7.015 

12.627 

(MHz) 

60 

2.355 

4.756 

4.874 

8.33 

0 

281 

378 

323 

454 

u!  *E= 

30 

328 

283 

375 

353 

60 

331 

261 

408 

346 

0 

77 

189 

107 

201 

ymF2 

30 

68 

105 

95 

177 

(km) 

60 

69 

81 

110 

137 
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Table  2-5.  E-  and  F-region  data  for  June,  40  degrees  longitude. 


[ Sunspot  Number 

10 

130 

Parameter 

Latitude 

Local  Time 

Local  Time 

0 

12 

12 

0 

0.222 

3.586 

EM 

4.458 

30 

0.199 

3.627 

m 

4.337 

60 

0.768 

3.071 

1.106 

3.653 

0 

3.661 

7.013 

8.069 

10.75 

30 

4.538 

6.517 

8.412 

9.66 

60 

3.421 

4.597 

5.191 

6.141 

As  in  the  current  model,  only  representative  day  (noon)  and 
representative  night  (midnight)  profiles  are  modeled.  Diurnal  model- 
ing greatly  increases  the  data  storage  requirements  for  the  E-  and 
F-regions  and  requires  detailed  energy  deposition  and  chemistry  model 
ing  in  the  D-region  where  most  of  the  variation  occurs  near  sunrise 
and  sunset. 

E-  and  F- Reg  ion  Model 

Use  of  the  parabolas  described  in  Equation  2-7  results  in 
discontinuities  in  electron  density  at  h^  + y . While  generally  not 
significant  for  the  intended  application,  it  is  desirable  to  prevent 
discontinuities  and  a profile  fitting  procedure  described  by  Booker 
(Reference  2-10)  is  used  to  model  the  E-  and  F-regions. 

In  order  to  prevent  discontinuities  in  the  electron  density 
or  any  of  its  derivatives,  Booker  suggests  fitting  the  electron  den- 
sity profile  with  exponential  fits  between  selected  altitudes  (transi 
tion  altitudes)  and  then  using  smoothing  functions  at  the  transition 
altitudes  which  smooth  the  electron  density  over  a prescribed  scale 
length.  The  relation  for  electron  density  presented  by  Booker  is 

log  Nc(z)  = log  Nq  + AQ1(z  - zq)  + 

1 (A  , - A , ) + 

n,n+l  n-l,n' 

n=l 

{f(,  - zn,Bn)  - f(zo  - Zn,Bn)}  (2-10) 

altitude 

transition  altitudes  used  to  describe  electron 
density  profile 

slope  of  log  Nq(z)  versus  z for  z^  < z < zn+j 

reciprocal  of  smoothing  scale  used  at  transition 
altitudes 


where 


z = h = 


n,n+l 


B 
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!B’*£n{l  + exp(Bz)}  zB  < 100 

z zB  > 100 

The  above  formulation  is  used  in  the  WEPH  code  model  by 
choosing  scaling  parameters  from  the  E-  and  F-region  parabola  param- 
eters given  in  Tables  2-4,  2-5,  and  2-6.  Figure  2-2  shows  an  example 
of  the  parabolas  used  to  fit  the  daytime  ionosphere.  Transition  alti- 
tudes used  in  determining  parameters  for  the  Booker  formulation  are 
shown  in  Figure  2-2  and  are  described  in  Table  2-7  along  with  the 
choice  for  the  reciprocal  of  B.  The  electron  density  slope  below 
Zj  is  found  by  setting  the  electron  density  at  h^  - y^  equal  to 

N^/1000.  The  electron  density  slope  above  Zy  is  found  by  computing  the 

the  electron  density  at  h^p^  + ^^2  ^Tom  Equation  2-8. 

Table  2-7.  Scaling  parameters  for 

use  in  Booker  formulation. 

Transition 


Altitude 

Description 

Smoothing  Scale 

Z1 

hmE  ' 0,8  yE 

0.5 

z2 

hmE 

0.1(yE) 

z3 

°'5(|lmE  + yE  + hmFl  “ yFl^ 

o.i (yE  + yF1 )/2 

Z4 

hmFl 

o.i (yF1) 

z5 

0 * 5( hmFl  + yFl  + hmF2  + yF2* 

o« i (y  pi  + yF2)/2 

z6 

hmF2 

o.i(yF2) 

z7 

hmF2  + yF2 

o.i (yF2) 
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Figure  2-2.  Parabolas  for  daytime  Ionosphere  and  selected 

transition  altitudes  for  use  in  Booker  formulation. 


After  the  electron  density  at  a specified  altitude  is  found, 
the  ion-pair  production  rate  is  determined  from  the  chemistry  model 
as  in  the  current  WEPH  code  ionosphere  model. 

D- Region  Model 

For  the  D-region  model  Equations  2-3  through  2-6  are  used  to 
define  the  electron  density  at  altitude  h^  and  h^  where 


h = h + — 
1 o a 


(2-11) 


w 


h_  = h - — 
2 o a 


(2-12) 


w 


Then  the  D-region  chemistry  model  is  used  to  determine  the  total 
effective  ion-pair  production  rates  (q^  and  q2)  required  to  produce 
the  specified  electron  densities  at  h^  and  h2.  In  computing  the  sun- 
spot number  the  following  relation  obtained  from  Reference  2-3  is 
used 

SSN  = - 406.35  + /93608  ♦ 1117  F 

where  F is  the  average  10.7  cm  solar  flux.  The  quantity  F is  obtained 
from  the  SAI  atmospheric  model  developed  for  ROSCOE  in  which  F is 
based  on  an  assumed  sinusodial  11-year  (4018  day)  variation  with  a 
maximum  value  of  250  and  a minimum  value  of  65. 

The  cosmic  ray  production  rates  at  hj  and  h2  are  found  from 
Equation  2-1,  and  ion-pair  production  rates  that  can  be  added  to  the 
cosmic  ray  production  rates  determined  from 


(\l  = q(hj)  - qcChj) 

q2  = q(h2)  - qc(h2) 


(2-14) 


(2-15) 


If  the  cosmic  ray  production  rate  is  larger  than  the  total  q required 
at  either  h^  or  h2>  the  cosmic  ray  production  rate  is  reduced  to  a 
value  slightly  less  than  the  total  production  rate  (q  = 0.9q). 
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Tor  daytime  conditions  an  effective  Lyman  a production  rate 
is  defined  in  terms  of  qj  and  q7  by 

qLa  = C[N0]exp(-6p)  (2- lb) 

where 


P(hj)  - p(h2) 


C 


ql 

[NO] : exp [- $p (hj ) ] 


[NO] j = nitric  oxide  density  at  hj 
[N0]?  = nitric  oxide  density  at 


The  Lyman  a and  cosmic  ray  production  rates  are  added  to  obtain  the 
total  production  rate.  The  above  procedure  results  in  a value  of  B 
for  a number  of  conditions  smaller  than  the  value  derived  for  Lyman 
radiation  in  Reference  2-7.  This  results  in  underestimating  the 
Lyman  a ionization  in  the  upper  D-region.  To  prevent  this  the  value 
of  B is  limited  to  be  equal  or  greater  than  0.05.  Then  to  maintain 
the  ionization  gradient  between  h^  and  h2  the  cosmic  ray  production 
rate  is  modified  so  that 


q2  = qLa(h2)  + qc(h2^  * (2'17) 

For  nighttime  conditions  an  exponential  ion-pair  production 
rate  is  determined  in  terms  of  qj  and  q2  from 

qE  = C exp[B(p) ] (2-18) 

where 
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1 


- 


' " exp [ - 3p (h 2 ) J 

Tlie  value  of  6 is  not  allowed  to  become  negative,  which  can  happen 
when  th  chemistry  and  electron  density  models  arc  not  consistent. 

Note  that  the  above  procedures  result  in  an  exponential  vari- 
ation in  production  rate  rather  than  electron  density  near  the  alti- 
tude h . However,  the  difference  in  the  electron  density  profile 
c 

near  h is  small  and  the  above  procedure  results  in  a better  fit  to 
the  electron  density  below  hQ  where  atmospheric  chemistry  can  affect 
the  gradient. 

Examples  of  Model  Output 

figures  2-3  and  2-4  show  comparisons  of  electron  density 
obtained  from  the  current  and  new  ionosphere  models  for  conditions 
similar  to  those  for  which  the  current  model  is  applicable. 

Figure  2-5  shows  the  D-region  results  from  the  new  model  in 
more  detail.  The  quantity  a is  the  reciprocal  of  the  electron  density 
scale  height.  Note  that  a is  not  constant  but  does  equal  specified 
in  Equation  2-6  near  hQ.  For  daytime  conditions  the  electron  density 
is  nearly  constant  between  hw  and  the  bottom  of  the  E-region.  The 
ion-pair  production  rate  and  electron  density  would  increase  above 
hw  if  the  NO  concentration  was  increased  (see  Section  1) . For  night- 
time conditions  the  altitude  h is  below  the  altitude  where  there  is 

o 

significant  change  in  the  electron  loss  rate  in  the  WEPH  code  chemistry 
model.  This  requires  an  increase  in  the  ion-pair  production  rate  with 
decreasing  altitude  and  results  in  a negative  value  for  8 in  Equation 
2-18.  Limiting  the  value  of  0 to  zero  results  in  a slightly  larger 
value  of  a near  h than  a . The  inconsistency  between  the  electron 
density  model  and  the  chemistry  model  could  be  resolved  if  more  de- 
tailed latitude,  season,  and  solar  cycle  variations  were  included  in 
the  chemistry  model.  The  limiting  value  of  0 also  results  in  a nearly 
constant  electron  density  between  hw  and  the  bottom  of  the  E-region 
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ELECTRON  DENSITY  (cm'J) 

e 2-3.  Comparison  of  current  and  new  ionosphere 
models  for  daytime  conditions. 


Figure  2-4.  Comparison  of  current  and  new  ionosphere 
models  for  nighttime  conditions. 


1 


Figure  2-5.  D-region  ionization,  30  degrees  latitude,  15  June  1961. 


for  nighttime  conditions.  Use  of  an  exponential  fit  between  and 
the  F.-rcgion  might  provide  a better  model. 

The  positive  and  negative  ion  densities  predicted  by  the 
model  also  depend  on  the  ion-pair  production  rate  and  the  atmospheric 
chemistry  model.  The  D-region  chemistry  model  used  in  the  WF.PH  code 
was  developed  for  moderate  to  large  production  rates  and  may  result 
in  significant  errors  when  used  to  predict  ambient  ion  concentrat- 
tions.  Also,  the  lack  of  latitude,  seasonal  and  solar  cycle  depen- 
dencies may  result  in  incorrect  variations  in  ion  densities  with  these 
parameters . 

Figures  2-6  and  2-7  show  day  and  night  profiles  obtained  with 
the  new  model  for  different  conditions.  The  variation  in  the  night- 
time F-region  electron  density  can  be  large  with  relatively  low  values 
occurring  at  midlatitudes  during  sunspot  minimum  (sec  Tables  2-4,  2-5, 
and  2-6) . 


New  Ion-Neutral  Collision  Frequency  Model 

As  previously  mentioned,  there  is  considerable  uncertainty  in 
ion-neutral  collision  frequency  predicitons  and  a simple  formulation 
in  terms  of  the  electron-neutral  collision  frequency  is  currently 
used  in  the  WFPH  code.  Recently,  the  U.S.  Army  Electronics  Command 
initiated  systematic  measurements  of  atmospheric  ion  mobility  (Refer- 
ence 2-11),  and  it  appears  appropriate  to  model  the  ion-neutral  colli- 
sion frequency  in  terms  of  ion  mobility. 

The  low  frequency  ion-neutral  collision  frequency  is  given  in 
terms  of  ion  mobility  in  Reference  2-12  as 


v. 

id 


9.6  x 10Un 

M u. 
r l 


-1 

sec 


(2-19) 


where 

n = neutral  particle  density,  cm 
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ELECTRON  DENSITY  (cut'3) 

Figure  2-6.  Comparisons  of  daytime  profiles  obtained 
with  new  ionosphere  model. 


OtUHBER  1%1 


Comparisons  of  nighttime  profiles  obtained 
with  new  ionosphere  model. 


M.M 

M = reduced  mass  = -1  -n-,r 

r M.  + M 

1 n 

M.  ,M  = ion  and  neutral  particle  mass  respectively,  AMU 

u.  = ion  mobility  V"*cm*sec  * 

1 1 

Note  that  the  reduced  mass  only  varies  by  a factor  of  two  as  the  ion 

mass  varies  from  M (~  30)  to  a value  much  larger  than  M . The  thermal 
n n 

ion  mobility  is  also  given  in  terms  of  the  reduced  mass  in  Reference 
2-12; 


yi 


1 

AT 


r 

Thus,  the  ion-neutral  collision  frequency  may 
tive  to  ion  mass. 


(2-20) 

be  relatively  insensi- 


For  an  assumed  ion  mass  of  60  AMU  and  vl  11  1//T  , liquation 


2-19  becomes 


v. 

tm 


9 1/2 

2.8  x 10'  nT  ' 


-1 

sec 


(2-21) 


Theoretical  calculation  for  ion  mobilities  in  air  given  in  Reference 

2-12  vary  from  about  9 x 10^V”*cm’ *sec"*  for  N0+  to  about  6 x lO1^ 

V ^cm  ^sec for  a hydrated  ion.  Atmospheric  measurement  for  positive 

and  negative  ion  mobilities  given  in  Reference  2-13  vary  from  about 
19  19  -1  -l  -l 

8 x 10  to  22  x 10  'V  cm  sec  (some  much  smaller  values  believed 
to  be  associated  with  aerosol  particles  are  also  reported  for  alti- 
tudes above  about  60  km).  The  above  range  in  ion  mobilities  results 

in  ion-neutral  collision  frequencies  that  vary  from  v /S  to  v /20 

n em  cm 

where  vem  is  the  high-frequency  electron-neutral  collision  frequency 

(v.  = \J  / 20  in  the  current  WEPH  code  model), 

lm  vcm 

In  performing  low  frequency  propagation  calculations  NOSC  has 
adapted  an  ion-neutval  collision  frequency  that  varies  exponentially 
from  lO^scc  * at  sea  level  to  8 x 10'\ec at  100  km  for  nighttime 
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conditions  and  is  about  half  these  values  for  daytime  conditions 

(Reference  2-14).  These  values  correspond  to  v /20  at  sea  level 

_ cm 

and  v /5  at  100  km.  The  smaller  values  of  v / in  the  upper  D- 
em  cm  em  ' 1 

region  appear  to  result  in  better  agreement  with  ambient  propagation 
data  than  the  larger  sea  level  values. 

The  following  expression  for  ion  mobility  is  used  in  Equa- 
tion 2-21  for  the  new  ion-neutral  collision  frequency  model 

!20  x 1019  h < 50  km 

20  x 1019  exp [-  0 . 277(h~50) ] 50  km  < h < 100  km 

5 x 1019  h > 1000  km  . (2-2 2) 

These  values  for  ion  mobility  are  reasonably  consistent  with  measure- 
ments reported  in  Reference  2-13  and  result  in  ion-neutral  collision 
frequencies  similar  to  those  used  at  NOSC  in  the  upper  D-region. 

COMPARISONS  WITH  PROPAGATION  DATA 
ELF  Propaaation  Data 

In  Reference  2-15,  78  Hz  propagation  measurements  between  Mt. 
Airy,  North  Carolina  and  Goose  Bay,  Labrador  (midpoint  ft  48°N  latitude, 
290°E  longitude)  are  presented  that  indicate  that  the  attenuation  rate 
varied  from  about  1 dB  per  megameter  at  night  to  about  1.4  dB  per  mega- 
meter  during  the  day.  Similar  values  were  presented  in  Reference  2-16 
for  an  eastwest  path  between  Hawaii  and  Utah.  Table  2-8  shows  attenu- 
ation rates  for  these  paths  obtained  with  several  ionosphere  and  colli- 
sion frequency  models. 

The  first  group  of  calculations  shown  in  Table  2-8  are  for 
the  north-south  path.  Figures  2-8  and  2-9  show  the  several  electron 
density  profiles  used.  The  nighttime  profile  labeled  WEDC0M  is  a modi- 
fication of  the  WEPH  profile  that  is  similar  to  the  profile  used  at 
NOSC  (Reference  2-14).  The  profiles  labeled  New  Model  (preliminary) 
were  obtained  with  a preliminary  version  of  the  new  ionosphere  model 
previously  described.  The  profiles  labeled  New  Model  with  NOSC 
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Table  2-8.  ELF  calculations. 


North  Carolina  to  Goose  Bay,  Labrador 
(Midpoints  48°  N Lat.,  290°  E Long.) 


Attenuation  Rate 


Ionosphere 

Model 

Collision 
Frequency  Model 

(dB/1000  km) 
Night  Day 

WEPH 

WEDCOM 

1.86 

1.64 

WEDCOM 

WEDCOM 

0.82 

1.64 

WEDCOM 

New 

1.0 

1.73 

New  (15  Dec  66) 

WEDCOM 

1.2 

1.22 

New  (15  Dec  66) 

WEPH 

1.46 

1.25 

New  (15  Dec  66) 

New 

1.42 

1.26 

New  (15  Dec  66) 
with  NOSC  parameters 

New 

0.78 

1.1 

Utah  to  Hawaii 
(Midpoint  a 30°  N Lat., 

220°  E Long.) 

New  (15  March  71) 

WEDCOM 

1.26 

- 

New  (15  March  61 ) 

WEDCOM 

1.19 

- 

New  (15  March  64) 

WEDCOM 

3.2 

_ 
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200 


Electron  density  profiles  used  for 
nighttime  ELF  calculations. 


I 70 


e 2-9.  Electron  density  profiles  used  for 
daytime  ELF  calculations. 


NEW  MODEL 


parameters  were  obtained  by  using  hw  = 88  km,  = 0.35  for  nighttime 
and  hw  = 72  km,  = 0.15  for  daytime  in  the  preliminary  version  of 
the  new  model.  The  profiles  labeled  New  Model  are  for  the  final 
version  of  the  new  ionosphere  model  and  are  shown  for  reference. 

They  were  not  used  for  low  frequency  calculations. 

Figure  2-10  shows  the  ion-neutral  collision  frequency  pro- 
files used.  The  profile  labeled  WEDCOM  is  similar  to  the  profile 
used  at  NOSC. 

The  daytime  results  are  relatively  insensitive  to  the  ion- 
neutral  collision  frequency  model,  but  are  sensitive  to  the  electron 
density  profile.  The  results  for  the  new  ionosphere  model  and  parti- 
cularly for  the  new  model  with  NOSC  parameters  are  more  consistent 
with  the  propagation  data  than  the  current  WEPH  and  WEDCOM  models. 

The  difference  appears  to  be  due  to  the  larger  gradient  in  electron 
density  between  80  and  100  km  in  the  new  model.  The  final  version  of 
the  new  model  has  a larger  gradient  than  the  preliminary  version 
shown  in  Figure  2-9  and  should  result  in  somewhat  lower  attenuation 
rates. 

For  nighttime  conditions,  both  the  electron  density  and  ion- 
neutral  collision  frequency  profiles  are  important.  The  WEDCOM  elec- 
tron density  profile  results  in  attenuation  rates  consistent  with  mea- 
sured values  with  either  the  WEDCOM  or  new  ion-neutral  collision  fre- 
quency models.  When  the  new  ionosphere  model  with  D-region  parameters 

(h  ,a  ) chosen  from  Equations  2-5  and  2-6  is  used,  the  attenuation 
w w 

rate  is  high  even  when  the  WEDCOM  ion-neutral  collision  frequency 
profile  is  used.  If  the  D-region  parameters  developed  from  propaga- 
tion data  by  NOSC  are  used,  the  nighttime  attenuation  rate  is  signi- 
ficantly reduced.  Calculations  are  also  shown  in  Table  2-8  for  a 
nighttime  path  between  Utah  and  Hawaii.  Except  for  sunspot  minimum 
(1964)  the  results  are  similar  to  those  for  the  northsouth  path. 

For  sunspot  minimum  conditions  near  30  degrees  latitude  the  E-region 
electron  density  is  predicted  to  be  quite  small  (see  Table  2-4  and 
Figure  2-7  and  the  resulting  ELF  attenuation  rate  is  large. 
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Clearly,  the  nighttime  attenuation  rates  are  critically  de- 
pendent on  the  details  of  the  electron  density  profile  in  the  D-  and 
E-regions.  Calculations  of  reflected  energy  similar  to  those  presented 
in  Reference  2-17  are  needed  to  study  the  relative  importance  of  the 
regions.  The  nighttime  electron  densities  obtained  with  the  new  model 
and  NOSC  parameters  are  smaller  than  would  be  predicted  from  the  de- 
cay of  daytime  electron  densities  with  the  WEPH  code  chemistry  model 
indicating  that  either  the  chemistry  model  is  in  error  (detachment 
rates  too  high  between  80  and  90  km),  or  the  NOSC  value  of  hw  is  too 
high.  The  nighttime  results  also  show  some  sensitivity  to  the  ion- 
neutral  collision  frequency  depending  on  the  ion  densities  in  the  D- 
region.  The  ion  densities  in  the  new  model  (particularly  the  new 
model  with  NOSC  0-region  parameters)  may  be  too  low  since  the  steady- 
state  value  of  ion  density  can  be  less  than  the  transient  decay  of 
the  daytime  ion  density.  Thus,  the  significance  of  the  ion-neutral 
collision  frequency  may  be  understated  if  the  actual  ion  densities 
are  larger  than  modeled. 

VLF  Propaqation  Data 

Figures  2-11  and  2-12  show  results  of  daytime  VLF  measurements 
for  propagation  between  Hawaii  and  San  Francisco  for  summer  conditions 
and  Hawaii  and  Utah  for  winter  conditions  presented  in  Reference  2-2. 
Also  shown  are  NOSC  calculations  for  exponential  electron  density  pro- 
files. As  discussed  in  Reference  2-2,  the  profile  defined  by  hw  = 70, 
aw  = 0.35  gives  the  best  results  for  summer  conditions  and  the  profile 
defined  by  hw  = 75,  a = 0.15  gives  the  best  results  for  winter 
conditions. 

WEDCOM  calculations  at  about  400  km  intervals  obtained  with 
the  new  ionosphere  (preliminary  version)  and  ion-neutral  collision 
frequency  models  are  shown  in  Figures  2-12  and  2-13  for  comparison 
with  the  measured  data.  The  results  for  the  24  kHz  propagation  dur- 
ing summer  are  in  reasonable  agreement  with  the  data  and  are  about 


AMPLITUDE,  dB  ABOVE  1 .V/m 
FOR  1-kW  RADIATED  POWER 
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Figure  2-11.  Propagation  from  San  Francisco  to  Hawaii, 

May  18-19,  1965  (NPM  24  kHz). 

midway  between  the  two  NOSC  curves  as  would  be  expected  from  the 

value  of  h and  a predicted  in  the  new  model, 
w w 

The  results  for  the  28  kHz  propagation  during  winter  show 
more  variation  from  the  data  and  from  the  NOSC  calculations  than  for 
the  previous  case.  Calculations  were  made  with  an  exponential  pro- 
file (hw  = 75,  aw  = 0.15)  to  check  the  propagation  calculations.  The 
results  are  essentially  the  same  as  the  NOSC  results  and  the  small 
differences  remaining  can  probably  be  attributed  to  differences  in 
magnetic  field,  propagation  angles  and  the  number  of  profiles  used 
to  define  propagation.  Calculations  were  also  made  with  the  new 
ionosphere  model  with  aw  = 0.15  and  hw  = 75  to  see  how  much  of  the 
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CURVES  FRO*  REFERENCE  2-2 


Figure  2-12.  Propagation  from  Utah  to  Hawaii 

(daytime,  February  1974),  (28  kHz) 


ELECTRON  DENSITY  (cm*3) 

Figure  2-13.  Daytime  electron  density  profiles 
used  for  VLF  calculations. 

difference  in  results  was  due  to  the  use  of  an  exponential  profile 
and  how  much  was  due  to  the  profile  parameters  hw  and  a^.  The 
several  electron  density  profiles  used  are  shown  in  Figure  2-13. 

For  the  case  shown  the  electron  density  profile  obtained  with  the 
preliminary  and  final  versions  of  the  new  model  are  essentially  the 
same.  It  would  appear  that  the  profile  parameters  recommended  by 
NOSC  provide  better  agreement  with  the  data  than  those  obtained  from 
Equations  2-5  and  2-6.  Note  that  the  field  strength  calculated  as 
the  vector  mode  sum  is  sensitive  to  relatively  small  changes  in  the 
electron  density  profile.  While  specific  profile  parameters  may  be 
chosen  emperically  for  ambient  conditions,  there  is  not  sufficient 
data  to  do  this  for  disturbed  conditions.  A more  useful  calculation 
for  disturbed  conditions  is  the  rms  mode  sum  which  is  not  as  sensi- 
tive to  the  details  of  electron  density  profile. 


Figure  2-14  shows  results  of  nighttime  VLF  measurements  pre- 
sented in  Reference  2-2  for  propagation  between  Hawaii  and  Ontario, 
California.  As  described  in  Reference  2-2,  nighttime  measurements 
show  variation  in  propagation  conditions  from  night  to  night  and 
throughout  a given  night.  WEDCOM  calculations  made  with  the  exponen- 
tial profile  do  not  show  the  same  interference  pattern  at  close-in 
distances  as  the  NOSC  calculations.  This  is  due  to  limiting  the  num- 
ber of  modes  considered  in  the  WEDCOM  propagation  model.  The  WEDCOM 
calculations  for  distances  beyond  3000  km  are  slightly  larger  than 
those  by  NOSC  and  would  be  even  larger  if  additional  modes  were  con- 
sidered. However,  the  attenuation  rate  is  sensitive  to  the  electron- 
neutral  collision  frequency  used  and  the  signal  strength  would  be  re- 
duced if  the  exact  collision  frequency  model  used  at  NOSC  was  used 
in  the  WEDCOM  calculations. 

The  calculations  for  the  new  ionosphere  model  with  h and  a 

w w 

determined  from  Equations  2-5  and  2-6  and  set  equal  to  NOSC  values 
result  in  larger  signal  strengths  than  the  exponential  model.  This 
is  due  to  the  larger  electron  density  gradient  near  the  reflection 
altitude  (see  Figure  2-15).  The  electron  density  gradient  is  larger 
for  the  new  model,  both  because  aw  is  larger  and  because  of  the  effect 
of  the  WF.PH  code  chemistry  model  as  previously  described.  When  the 
NOSC  values  for  hw  and  aw  are  used  in  the  new  model,  the  total  pro- 
duction rates  are  less  than  the  nominal  cosmic  ray  values,  and  the 
method  of  reducing  the  cosmic  ray  ion  production  rate  and  determining 
the  total  production  rate  results  in  a larger  electron  density  gradi- 
ent. The  final  version  of  the  new  model  results  in  an  electron  den- 
sity profile  similar  to  that  for  the  preliminary  model  when  the  NOSC 
parameters  are  used. 

Nighttime  calculations  at  the  upper  end  of  the  VLF  band  are 
the  most  difficult  to  do  accurately  and  are  affected  by  relatively 
small  changes  in  the  electron  density  gradient.  As  for  daytime  cal- 
culations, the  vector  mode  sum  is  particularly  sensitive  to  changes 
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Figure  2-15.  Nighttime  electron  density  profile  used 
for  VLF  calculation. 

in  the  electron  density  profile  and  requires  detailed  matching  to 
predict  measured  values  (and  consideration  of  a relatively  large  num- 
ber of  modes).  The  rms  mode  sum  is  less  sensitive  and,  while  more 
variable  with  electron  density  profile  parameters  than  the  daytime 
results,  gives  estimates  of  expected  signal  strength  that  are  within 
the  variations  occurring  from  night  to  night. 

HF  Propagation  Data 

An  ambient  HF  absorption  relation  is  used  in  the  ITS  predic- 
tion code  that  is  based  on  propagation  data  (Reference  2-18).  When 
converted  to  one-way  vertical  absorption  for  12:00  hours  local  time 
the  relation  becomes 
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(2-23) 


A _ 143(1  + 0.0087|e|)W 
10.2  + (f  + 1 . 4) 2 

• ((1  + 0.005  S)  cosmx] 


where 


W 


1.807  - 0.0269  | 1 9 | - 60 1 for  winter  and  1 0 1 5 30 
1 for  other  conditions 


m = maximum  (1.3,  2.25  - f* . 032 1 9 1 ) 
f = frequency  (MHz) 

S = smoothed  1 3-month  sunspot  number 
x = solar  zenith  angle  at  noon 
0 = geographic  latitude  (degrees) 


The  ITS  prediction  program  is  being  revised  and  a new  absorp- 
tion relation  (Reference  2-19)  developed.  A preliminary  version  of 
the  relation  (converted  to  one-way  vertical,  noon  time  absorption)  is 


A = 


338  I 

10.2  ♦ (f  + 1.4) 2 


(2-24)  „ 


where 


I = - 0.04  ♦ exp (-2. 9 37  + 0.8445  fQE)  . 

In  this  expression  the  relation  for  absorption  index  I is  in  terms  of 
fgp,  which  includes  the  variation  in  latitude  and  solar  activity. 

Table  2-9  shows  comparisons  of  the  predicted  and  calculated  one-way 
vertical  absorption  at  15  Mlz  for  different  latitudes,  seasons,  and 
sunspot  numbers. 

The  absorption  results  obtained  with  the  new  model  (final  ver- 
sion) are  generally  between  those  obtained  with  Equations  2-23  and 
2-24.  The  summer  predictions  tend  to  be  low  except  for  60  degrees 
latitude  where  the  prediction  is  high.  Changes  in  the  NO  concentration 
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in  the  upper  D- region  would  affect  the  results  as  well  as  changes  in 
the  scaling  for  hw  and  a . Since  the  electron  density  profile  in  the 
new  model  is  dependent  on  both  the  E-rcgion  critical  frequency  and 
the  l)-rcgion  parameters  h^  and  a , the  variation  with  latitude  and 
season  is  different  than  the  new  ITS  model,  which  is  only  dependent 
on  the  E-region  critical  frequency. 

SUMMARY 

The  structure  of  the  new  ionosphere  model  provides  for  lati- 
tude, season,  and  solar  cycle  variations.  A better  data  base  for  D- 
region  electron  density  profiles  is  needed  to  improve  calculations  of 

h and  a . Partial  reflection  measurements  being  made  on  a routine 
w w b 

basis  at  several  latitudes  by  the  U.S.  Army  Electronics  Command  may 
be  useful  in  providing  the  data  base. 

Values  for  hw  and  obtained  by  NOSC  from  propagation  mea- 
surements can  be  used  in  the  model  instead  of  values  obtained  from 
fits  to  electron  density  data  (Equations  2-5  and  2-6).  This  would 
improve  comparisons  with  measured  VLF  data  presented  by  NOSC,  parti- 
cularly for  daytime  conditions.  Additional  scaling  parameters  to 
aid  in  modeling  the  electron  density  between  hw  and  the  bottom  of  the 
E-region  are  needed.  The  D-region  predictions  would  also  be  improved 
by  having  the  atmospheric  chemistry  model  dependent  on  latitude,  sea- 
son, and  solar  conditions,  llse  of  the  current  WE  PH  code  model  may 
result  in  errors  in  predicting  ion  densities  for  both  day  and  night 
conditions  and  the  electron  density  gradient  near  80  km  for  nighttime 
conditions.  Use  of  steady-state  relation  may  also  underestimate  night- 
time electron  and  ion  densities. 

Further  study  of  ELF  propagation  data  is  needed  to  determine 
the  critical  aspects  of  the  F.- region  profile  and  whether  the  predicted 
profiles  are  adequate  for  propagation  calculations.  For  both  day  and 
night  conditions  the  VLF  vector  mode  sum  is  critically  dependent  on 
the  electron  density  profile  and  the  rms  mode  sum  is  the  more  useful 
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descriptor  for  weapon  effects  codes.  Benchmark  data  for  comparison 
of  HF  absorption  predictions  are  needed  to  determine  whether  the 
model  is  adequate  for  propagation  calculations. 

The  major  change  to  the  ion-neutral  collision  frequency 
model  was  to  parameterize  the  model  in  terms  of  ion  mobility.  While 
ambient  propagation  measurements  show  some  sensitivity  to  the  ion- 
neutral  collision  frequency  model,  uncertainties  in  electron  and  ion 
density  profiles  prevent  further  refinement  of  the  ion  mobility 
expression. 
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SECTION  3 
HEAVE  MODEL 


INTRODUCTION 

The  fireball  motion  and  atmospheric  heave  models  are  uncoupled 
in  the  current  WF.PH  code  and  the  heave  model  can  underestimate  atmos- 
pheric motion  near  the  burst  point.  Thus,  it  is  possible  to  have  a 
highly  heated  and  highly  ionized  region  below  the  fireball  that  should 
have  risen  with  the  fireball.  Because  of  fircbal 1-likc  properties, 
this  region  has  been  called  a "ghost  fireball".  In  the  RANC  code,  the 
ghost  fireball  problem  was  minimized  by  enlarging  thT*  initial  fire- 
ball size  to  include  air  that  could  remain  highly  ionized.  However, 
when  the  ROSCOF.  energy  deposition  models  are  used  in  conjunction  with 
RANC  fireball  models,  the  initial  fireball  size  is  too  small  to  pre- 
vent a ghost  fireball  for  certain  burst  conditions. 

The  atmospheric  heave  model  in  the  WF.PH  VI  code  was  developed 
for  the  WOF.  code  to  reduce  computation  time  relative  to  the  previous 
heave  model  developed  for  thj  RANC  code.  The  WOE  code  model  also 
allows  for  "tying"  the  fireball  and  upwelling  atmosphere  together  in 
terms  of  a fireball  wake  model  (Reference  3-1).  The  fireball  wake 
model  forces  the  ghost  fireball  into  the  same  region  occupied  by  the 
real  fireball. 


Fireball  Wake  Model 

The  fireball  wake  is  described  by  a right  circular  cylinder 
with  radius 


/ Ru  + RD  / ) 

Rw  = minimum  ^10 ^ , /1.5  R^  f 


(3-1) 
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where 


R^.Rp  = upper  and  downward  initial  fireball  radii, 
respectively 

Rw  = fireball  magnetic  radius 
M 

The  top  of  the  wake  column  is  determined  as  the  maximum  of  either  the 
altitude  a sonic  wave  would  reach  (assumed  velocity  is  0.5  km  s ) or 
one-half  of  the  upward  fireball  radius  beyond  the  fireball;  this 
modeling  produces  results  somewhat  similar  to  the  bow  wave  that  forms 
around  a blunt  object  at  subsonic  velocities  in  much  denser  air. 

The  wake  model  is  initialized  by  obtaining  the  Lagrangian 
coordinates  of  the  top  and  a point  near  the  bottom  of  the  initial  fire- 
ball. The  Lagrangian  coordinate  of  a point  is  the  mass  inventory  of 
the  air  column  above  that  point: 


oo 


After  some  variations,  a point  one  scale  height  below  the  bottom  of 
the  initial  fireball  was  chosen  as  the  bottom  boundary  of  the  rising 
fireball  for  the  WBPH  code  model.  This  encompasses  air  that  can  be 
highly  heated,  dissociated,  and  ionized  by  the  initial  energy 
deposition. 

A one-dimensional  model  allowing  only  vertical  flow  maintains 
the  Lagrangian  coordinate  of  every  air  parcel  as  an  invariant  property. 
By  adjusting  the  computed  Lagrangian  coordinates  in  the  burst  axis 
air  column  to  fit  the  predicted  fireball  position  and  its  initially 
obtained  coordinates,  the  developing  fireball  "draws"  external  air  up 
with  it  as  it  rises.  Figures  3-1  and  3-2  schematically  indicate  the 
initial  definition  of  the  fireball’s  Lagrangian  coordinates  and  their 
subsequent  use  to  adjust  the  heave  model  predictions  in  the  vicinity 
of  the  fireball. 
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FIREBALL  WAKE  CYLINDER 


A weighting  function  controls  the  adjustment  of  Lagrangian 
coordinates  horizontally  within  the  wake  cylinder.  It  is  a function 
of  the  ratio  of  horizontal  range  from  the  cylinder  axis  and  the  cylin- 
der's radius: 

K - JT  (5-3) 

1 + exp[16(r/Rw) 

where  r is  the  horizontal  range.  Vertically,  a linear  adjustment  is 
made  for  points  above  the  top  of  the  fireball  or  below  the  bottom 
boundary  (taken  as  one  scale  height  below  the  fireball  bottom).  For 
points  between  the  bottom  boundary  and  the  top  of  the  fireball  the 
Lagrangian  coordinate  value  is  scaled  to  match  the  boundary  values 
and  to  approximate  air  motion  within  the  fireball  region.  However, 
the  model  is  not  intended  for  detailed  modeling  of  the  fireball  air 
properties  and  has  been  adjusted  to  underestimate  electron  density 
within  the  fireball  region.  The  fireball  overlay  model  provides 
estimates  of  fireball  electron  density. 

Up  to  10  fireball  wakes  may  be  superimposed.  When  overlaps 
occur,  bursts  are  treated  in  their  detonation  order.  In  a multiburst 
scenario,  Lagrangian  coordinates  arc  determined  by  the  composite 
effect  of  atmospheric  heave  and  existing  fireball  wakes.  The  coordi- 
nates of  each  successive  fireball  arc  determined  in  the  wake  environ- 
ment of  all  previous  fireballs. 

Table  3-1  shows  a comparison  of  electron  density  versus  alti- 
tude along  the  vertical  through  the  burst  point  240  seconds  after  a 
multimegaton  burst  detonated  in  the  E-rcgion  with  and  without  the  wake 
model.  The  bottom  of  the  initial  fireball  (on  the  burst  axis)  is  at 
104  km  and  the  bottom  of  the  fireball  at  240  seconds  is  at  486  km. 

Also  shown  in  Table  3-1  is  the  starting  altitudes  for  air  reaching  a 
given  altitude  at  240  seconds.  When  the  wake  model  is  not  used  the 
air  starts  at  points  within  the  initial  fireball  and  even  though  the 
energy  deposition  is  modified  to  not  be  larger  than  that  at  the  fire- 
ball edge,  the  air  is  dissociated  and  remains  highly  ionized.  When 
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Tabl«y3-1.  Comparison  of  electron  density  and  initial  air 

altitude  at  240  seconds  after  burst  on  the  burst 
axis  with  and  without  wake  model. 

No  Wake  Wake 


Altitude 

(km) 

Electron  Density 
(cm- 3) 

Initial 

Altitude 

(km) 

Electron  Density 
(cm"3) 

Initial 

Altitude 

(km) 

120 

6.9(10) 

107.8 

3.8(4) 

93.2 

140 

3.4(10) 

110. 2 

2.6(4) 

93.5 

160 

2.6(10) 

112.3 

2.2(4) 

93.7 

180 

2.0(10) 

114.4 

1.8(4) 

93.8 

200 

1.6(10) 

116.3 

1.2(3) 

93.9 

220 

1.3(10) 

118.3 

9.3(2) 

93.9 

240 

1.2(10) 

120.0 

7.4(2) 

94 

260 

1.1(10) 

122.2 

6.0(2) 

94 

280 

1.1(10) 

124.1 

6.0(2) 

94 

300 

9.7(9) 

126.2 

6.0(2) 

94 

the  wake  model  is  used  the  air  starts  far  enough  below  the  fireball 
that  ionization  is  quickly  lost  by  dissociative  recombination. 

A late- time  ghost  fireball  can  also  be  caused  by  procedures 
used  in  modeling  air  motion  and  changes  in  air  density  for  times  after 
burst  when  air  has  returned  to  its  initial  altitude.  In  the  current 
WE: PI i code  model  early-time  motion  is  estimated  by  computing  the  air 
location  and  density  at  100-sccond  intervals  starting  at  300  seconds 
after  the  last  burst  and  continuing  until  100  seconds  after  the  first 
burst.  While  generally  providing  good  estimates  of  the  air  motion, 
the  procedure  docs  miss  the  apogee  time  sufficiently  in  some  case  to 
significantly  underestimate  the  amount  of  volume  expansion  during 
rise.  A new  procedure  in  which  the  apogee  times  between  bursts  and 
after  the  last  burst  are  estimated  has  been  prepared  that  provides 
much  better  estimates  of  apogee  times  and  thus  the  maximum  effects  of 
volume  expansion. 


After  reaching  apogee  the  air  falls  and  the  air  density  in- 
creases. In  the  WEI’ll  code  model  the  increase  in  air  density  is  inter- 
preted as  air  entrainment.  This  procedure  is  intended  to  return  air 
to  preburst  conditions  when  the  heave  motion  has  ended.  In  the  cur- 
rent model  the  properties  of  the  entrained  air  are  taken  as  those  of 
the  ambient  air  at  the  altitude  of  the  air  parcel  as  it  falls.  This 
procedure  tends  to  underestimate  air  entrainment  near  the  final  alti- 
tude and  in  some  cases  the  air  temperature  remains  above  ambient  after 
reaching  the  final  altitude.  In  the  revised  model  the  properties  of 
the  entrained  ail  are  always  taken  as  those  of  ambient  air  at  the 
final  altitude. 

The  heave  model  was  developed  to  provide  estimates  of  changes 
in  air  density  and  interpretation  of  the  model  in  terms  of  particle 
motion  is  fictitious,  particularly  after  the  air  parcels  reach  apogee. 
The  current  modeling  including  the  use  of  entrainment  is  a method  for 
using  the  heave  model  in  a manner  that  is  expected  to  provide  reason- 
able estimates  for  a few  hundred  seconds  (until  apogee  times)  and  that 
will  return  the  atmosphere  to  ambient  conditions  at  late  times.  How- 
ever, since  the  motion  model  is  independent  of  the  air  temperature, 
there  can  still  be  cases  where  air  temperatures  remain  above  ambient 
at  late  times.  These  cases  occur  when  the  air  motion  due  to  X-rays 
and  heavy  particle  energy  deposition  is  small  and  heating  due  to  UV 
energy  deposition  raises  the  gas  temperature.  It  can  also  occur  when 
energy  in  metastahlc  states  is  returned  to  kinetic  energy  and  heats 
the  gas.  Generally,  the  increase  in  air  temperature  above  ambient  is 
relatively  small  for  these  cases. 
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SECTION  4 
CHEMISTRY  MODELS 


INTRODUCTION 

Outside  the  fireball,  separate  atmospheric  chemistry  models 
are  used  in  the  KEPH  VI  code  for  the  D-region  (altitudes  below  100 
km)  and  the  E-  and  F-regions  (altitudes  above  100  km).  Because  of 
simplifications  used  in  the  models,  there  can  be  discontinuities  in 
ionization  near  the  boundary  between  the  two  models.  While  this  is 
allowed  for  in  some  calculations  (eg,  refraction),  it  can  result  in 
potential  errors  in  calculating  propagation  effects.  The  following 
describes  chemistry  model  improvements  to  minimize  ionization  dis- 
continuities and  to  incorporate  new  chemistry  modeling  recently  pre- 
pared for  ROSCOE. 

D- REGION  CHEMISTRY  MODEL 

The  D-region  chemistry  model  used  in  the  WEPH  VI  code  was 
first  developed  for  ROSCOE  (Reference  4-1).  For  daytime  conditions 
the  effect  of  burst  produced  odd  nitrogen  on  the  ambient  values  of 
0,  03,  02(*A),  NO,  and  NO^  is  estimated  from  the  total  odd  nitrogen 
produced  and  effective  photodissociation  rates  determined  from  the 
preburst  ambient  species  concentrations.  The  ROSCOE  chemistry  models 
have  recently  been  updated  to  include  species  required  in  optical  pre- 
dictions (Reference  4-2) . In  performing  this  work  it  was  found  that 
including  the  odd  hydrogen  species  H,  OH,  and  H02  in  the  model  could 
significantly  modify  the  effect  of  burst-produced  odd  nitrogen  on  the 
late-time,  ambient  daytime  atmosphere. 

In  order  to  include  the  new  species  in  the  WEPH  code  it  is 
first  necessary  to  add  the  species  to  the  ambient  atmosphere  model. 
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The  ambient  model  prepared  for  ROSCOF.  by  SA1  (Reference  4-3)  has  been 
used  to  modify  the  WITH  code  ambient  model.  Table  4-1  shows  the  day 
and  night  ambient  values  for  the  odd  hydrogen  species  used  in  the 
model.  The  species  concentrations  arc  not  functions  of  location,  sea- 
son, or  solar  conditions.  Next,  the  relations  determining  the  effec- 
tive photodissociation  rates  are  modified  to  include  the  new  species 
and  reactions  involving  the  new  species  are  added  to  those  determining 
the  effect  of  burst-produced  odd  nitrogen.  In  determining  the  new 
ambient  species  concentrations  and  the  transient  decay  from  burst- 
produced  conditions  the  formulations  given  in  Reference  4-1  have  been 
modified  to  simplify  the  calculations  and  to  insure  convergence  of 
the  solutions  (Reference  4-2). 

Figure  4-1  shows  the  change  in  the  daytime  ambient  odd  oxygen 
species  (0  + 0,)  versus  the  amount  of  odd  nitrogen  (NO  + NO,,)  . The 
subscript  o refers  to  preburst  ambient  conditions.  The  two  curves 
shown  in  Figure  4-1  labeled  WITH  code  are  calculations  for  the  WITH 
code  neutral  species  model  (developed  by  SAI  for  ROSCOF)  and  show  the 
effect  of  including  and  neglecting  the  odd  hydrogen  species.  While 
including  the  odd  hydrogen  species  has  considerable  effect  on  the 
reduction  of  the  oxygen  species,  the  effect  is  less  than  expected 
based  on  benchmark  calculations  prepared  by  Schiebc  for  development 
of  the  new  chemistry  models  for  ROSCOF;.  The  two  curves  in  Figure  4-1 
labeled  ROSCOF  Benchmark  are  for  the  ambient  neutral  species  used  by 
Scheibe  in  preparing  the  benchmark  calculations.  Note  that  the  reduc- 
tion in  odd  oxygen  is  much  less  for  Scheibe's  neutral  species  model 
even  when  the  odd  hydrogen  species  arc  neglected.  In  order  to  de- 
termine the  cause  of  the  difference,  calculations  were  made  for  a 

10  - 3 

weapon-produced  odd-nitrogen  concentration  of  4 x 10  cm  and 
changes  were  made  in  the  preburst  distribution  of  oxygen  and  nitro- 
gen species.  Figure  4-2  shows  calculation  as  a function  of  the  ratio 

of  [0]  to  £ where  £ = [0]  + [0_]  was  held  constant  at  2 x 1010 

^ Jo  o o 1 'o  1 3Jo 

cm  . Figure  4-3  shows  similar  results  as  a function  of  the  ratio  of 
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Table  4-1. 

Ambient 

species 

concentrations  (cm  ) 

for  H, 

OH,  and  HO2  used  in 

ROSCOE. 

Altitude 

H 

OH 

HOo 

(km) 

Day 

Night 

Day 

Night 

Day 

Night 

0 

7 . 0( -3 ) 

0.0 

1.0(6) 

1.7(2) 

1.0(5) 

4.9(1) 

5 

7.6(-3) 

0.0 

1.0(6) 

1.8(2) 

7.5(5) 

4.2(2) 

10 

1 . 0 ( -3 ) 

0.0 

1.05(6) 

2.1(2) 

2.4(6) 

1.6(3) 

15 

1.6(-2) 

0.0 

1.15(6) 

2.7(2) 

6.9(6) 

5.9(3) 

20 

5. 2(-2) 

0.0 

1.5(6) 

4.2(2) 

1.15(7) 

1.4(4) 

25 

3.2(-l ) 

0.0 

2.3(6) 

8.1(2) 

1.5(7) 

2.7(4) 

30 

2.9(0) 

0.0 

4.0(6) 

2.0(3) 

1.6(7) 

4.7(4) 

35 

1.0(2) 

0.0 

6.8(6) 

8.0(3) 

1.5(7) 

8.3(4) 

40 

4.0(4) 

0.0 

1.05(7) 

5.7(4) 

1.2(7) 

1.3(5) 

45 

1.0(5) 

0.0 

1.1(7) 

2.9(5) 

9.1(6) 

2.4(5) 

50 

2.4(5) 

0.0 

9.5(6) 

1.2(6) 

6.6(6) 

4.6(5) 

55 

5.1(5) 

0.0 

7.2(6) 

4.4(6) 

4.2(6) 

6.9(5) 

60 

1.0(6) 

0.0 

5.3(6) 

6.5(6) 

2.2(6) 

7.3(5) 

65 

1.8(6) 

0.0 

3.7(6) 

5.9(6) 

7.9(5) 

4.6(5) 

70 

4.9(6) 

0.0 

2.5(6) 

4.5(6) 

4.2(6) 

3.5(6) 

75 

1.25(7) 

5.0(2) 

1.6(6) 

3.2(6) 

1.2(7) 

1.2(7) 

80 

2.5(7) 

1.0(8) 

7.0(5) 

1.6(6) 

9.2(6) 

9.2(6) 

85 

8.6(7) 

8.6(7) 

7.0(4) 

1.7(5) 

5.7(4) 

5.7(4) 

90 

7.4(7) 

7.4(7) 

6.3(3) 

1.7(4) 

5.7(3) 

5.7(3) 

95 

5.0(7) 

5.0(7) 

5.7(2) 

1.7(3) 

4.9(2) 

4.9(2) 

100 

3.0(7) 

3.0(7) 

6.7(1) 

2.2(2) 

7.4(1) 

7.4(1) 
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Changes  in  the  steady-state  odd  oxygen  concentration 
at  60  km  as  a function  of  the  odd  nitrogen  concentration 


Figure  4-2.  Changes  in  the  steady-state  odd  oxygen  concentration  as 
a function  of  [0]  /£  , 6 = 4 x 101*  cm-3. 


Figure  4-3.  Changes  in  the  steady-state  odd  oxygen 
cohcentration  at  60  km  as  a function 
of  [N02]o/6o,  6 = 4 x 1010  cm"3. 
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the  ratio  of  fNCLl  to  6 where  fi  = fNO-1  + [N01  was  held  constant 
at  1.3  x lfT  cm* ^ . The  calculations  show  that  the  reduction  in  odd 
oxygen  caused  by  hurst-produced  odd  nitrogen  is  very  sensitive  to  the 
preburst  ratio  of  fN0o]o  to  S . This  is  due  to  the  computation  of 
effective  photodissociation  rates  from  the  preburst  ambient  species. 

In  the  Vs  I PH  code  model  the  ambient  NO,,  density  at  60  Km  is  large 
enough  to  cause  the  effective  photodissociation  of  N0?  to  become  nega- 
tive. further  detailed  studies  with  mill  t i spcci  cs  chemistry  codes  arc 
needed  to  determine  the  critical  parameters  affecting  the  reduction 
in  oxygen  species  and  the  allowable  range  in  these  parameters.  A 
better  way  of  determining  the  preburst  ambient  species  is  also  needed 
that  is  consistent  with  the  chemistry  models  used.  Above  about  50  km 
the  daytime  concentration  of  NO.,  is  determined  by  photodissociation 
and  the  ambient  daytime  value  should  be  limited  by  the  photodissocia- 
tion rate,  figure  4-4  shows  the  current  WfPH  code  daytime  ambient 
values  for  NO.,  and  modified  values  determined  by  computing  the  NO.,  con 
centration  from  the  photodissociation  rate  used  by  Schcibc  when  the 
effective  photodissociation  rate  determined  from  the  species  concen- 
trations becomes  less  than  that  used  by  Scheibe.  Use  of  the  modified 
ambient  NO.,  concentration  results  in  a much  smaller  reduction  in  odd 
oxv  gen  concentration  at  6(1  km  for  a given  burst-produced  odd  nitrogen 
concentration  as  shown  in  figure  4-3.  figure  4-5  shows  the  reduction 

in  daytime  odd  oxygen  as  a function  of  altitude  due  to  the  odd  ni tro- 
ll -3 

gen  produced  by  an  ionization  impulse  of  10  ion  pairs  cm  ‘ . Since 
most  of  the  daytime  detachment  rate  is  due  to  oxygen  species,  the 
amount  of  reduction,  particularly  above  about  55  km,  can  be  signifi- 
cant in  propagation  predictions. 

In  the  current  WfPH  code  the  nitrogen  vibrational  temperatures 
and  electron  kinetic  temperature  are  not  modified  by  energy  deposition 
in  the  II- region  chemistry  model.  This  causes  a discontinuity  in  the 
reaction  rate  coefficients  (primarily  the  dissociative  recombination 
coefficient)  and  thus  the  electron  density  near  100  km.  In  addition 
to  the  odd  hydrogen  species,  the  new  D- region  chemistry  model  for 
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Figure  4-5.  Reduction  In  odd  oxygen  as  a function  of 
altitude,  Nq  = 10n  ion  pairs  cm"3. 


ROSCOE  includes  a prediction  for  the  concentration  of  the  first  vibra- 
tional level  of  molecular  nitrogen.  If  a Boltzmann  distribution  for 
the  number  density  in  the  excited  state  is  assumed,  a nitrogen  vibra- 
tional temperature  can  be  computed  from 

Tv  = ev  (4-1) 


where 


B = 0.5  - (0.25  - N2(v=l)/N2r'‘ 

N2(v=l)  = concentration  of  molecular  nitrogen  in 
first  vibrational  level 


= total  concentration  of  molecular  nitrogen 
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If  it  is  also  assumed  that  the  nitrogen  vibrational  and  elec- 
tron kinetic  temperatures  are  equilibrated,  then 

Tc  = Tx  (4-2) 

where  Te  is  the  electron  temperature.  A modified  dissociative  recom- 
bination coefficient  can  be  computed  from 

ad  = “dT1  C4'3) 

e 

where 

a.  = dissociative  recombination  coefficient  computed 
d 

assuming  Te  = and  Tg  equals  the  gas  kinetic 
temperature  . 

Figure  4-6  shows  T versus  time  for  parametric  D-region  altitudes 

* 11  -3 

following  an  energy  deposition  impulse  producing  10  ion  pairs  cm  . 

A number  of  the  D-region  rates  have  been  modified  since  prepa- 
ration of  the  chemistry  model  for  ROSCOE.  Appendix  A shows  rates 
given  in  Reference  4-1  (used  in  the  current  WEPH  code)  and  revised 
rates  suggested  by  Niles  (Reference  4-4)  and  Scheibe  (Reference  4-5). 

E-  AND  F-REGION  CHEMISTRY  MODEL 

The  E-  and  F- region  chemistry  model  used  in  the  WEPH  code 
was  developed  for  ROSCOE  and  includes  modeling  of  an  excitation  tem- 
perature that  describes  equilibrium  Boltzmann  distributions  of  free 
electron  kinetic  energy,  oxygen  electronic  states,  and  nitrogen  vibra- 
tional states.  The  excitation  energy  is  modified  due  to  the  changing 
species  concentrations,  by  elastic  collisions  between  electrons  and 
heavy  particles,  quenching  of  the  O^D)  by  N2,  quenching  of  the  N2 

vibrational  states  by  CO.  and  heating  or  cooling  of  electrons  during 

^4 

collisional  excitation  of  N(  S)  of  deexcitation  of  N(4’D).  In  order 
to  improve  the  interface  between  the  E-  and  F-region  chemistry  model 
and  the  D-region  chemistry  model,  quenching  of  N2  vibrational  states 
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Figure  4-6.  Nitrogen  vibrational  temperature  following 
ionization  impulse  of  1011  ion  pairs  cm-3. 


by  0 and  09  has  been  added  and  the  modeling  of  collisional  excitation 
and  dccxcitation  of  nitrogen  atoms  has  been  modified.  Also,  the  forma- 
tion of  N(S)  and  N("D)  by  X-rays  has  been  modified  and  reactions  affec- 
ting the  formation  and  loss  of  NO  from  nitrogen  atoms  have  been  revised 

2 

to  account  for  the  rapid  formation  of  NO  by  N(  D) . 

In  the  h-  and  F- region  model  the  initial  species  concentrations 
are  determined  from  relations  developed  by  SAI  for  X-rays,  heavy  parti- 
cles, and  UV  (Reference  4-6).  For  X-rays  the  change  in  species  due  to 

the  initial  dissociation  is  neglected.  In  the  O-region  model  the 

4 

effects  of  dissociation  arc  included  and  result  in  0.45  N(  S)  atoms 

2 7 

and  0.61  N(  0)  atoms  per  ion  pair.  Neglect  of  the  production  of  N(“D) 

4 

by  dissociation  can  significantly  affect  the  decay  of  N(  S)  near  100  km. 
The  relations  used  below  100  km  have  also  been  used  above  100  km  to 
account  for  quenching  of  N9  vibrational  energy.  Since  water  molecules 
arc  not  modeled  above  100  km,  reactions  involving  water  molecules  in 
Archer's  model  were  eliminated  (Reference  4-2).  Figure  4-7  shows  a 
comparison  of  the  nitrogen  vibrational  temperature  computed  at  100  km 
altitude  with  the  D-region  model  and  the  E-  and  F-region  model  follow- 
ing a particular  energy  deposition  impulse.  The  initial  decay  of  vi- 
brational temperature  is  nearly  the  same  for  both  models.  The  ambient 
gas  temperature  at  100  km  is  200°K.  When  the  0- region  model  is  used 
the  vibrational  temperature  decays  until  reaching  about  220°K.  This 
temperature  is  due  to  the  production  of  vibrational  energy  by  the  ambi- 
ent ion-pair  production  rate.  When  the  E-  and  F-rcgion  model  is  used 
the  vibrational  temperature  decays  until  reaching  about  400°K.  This 
temperature  is  the  gas  kinetic  temperature  which  has  been  raised  in 
the  E-  and  F-region  model  due  to  energy  deposition.  The  heave  model' 
does  not  predict  enough  motion  to  result  in  cooling  the  air  and  it 
remains  above  ambient. 

The  calculation  procedure  used  in  the  current  E-  and  F-region 
chemistry  model  is  to  first  treat  the  ion  chemistry  (reactions  1 
through  16  in  Appendix  B) , including  changes  in  ion  concentrations 
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Figure  4-7.  Comparison  of  nitrogen  vibrational  temperature  at 

interface  between  D-region  and  E-  and  F-region  models. 


and  the  associated  changes  in  neutral  species  concentrations.  Then 
the  neutral  chemistry  (reactions  19  through  25)  is  computed  using 
effective  mean  values  of  the  species  determined  from  the  result  of 
the  ion  chemistry.  The  mean  values  are  computed  assuming  an  expo- 
nential change  in  species  concentration  with  time  during  the  time 
interval.  It  has  been  found  that  the  assumption  of  an  exponential 
change  can  overestimate  the  mean  value  of  electron  density  when  dis- 
sociative recombination  dominates  the  early  decay  of  ionization. 

Also,  the  use  of  mean  values  can  overestimate  reactions  between  species 
produced  during  the  time  interval. 

In  order  to  improve  the  neutral  chemistry  calculations,  parti- 
cularly near  the  boundary  between  the  D-region  and  E-  and  F-region 
models,  the  order  of  the  neutral  chemistry  calculations  and  the  approx- 
imations used  have  been  modified.  In  the  current  model  the  nitric 

oxide  chemistry  (reactions  5,  19  through  22,  and  25)  is  computed  first 

4 2 

and  then  the  equilibration  of  N(  S)  and  N(  D)  due  to  reactions  23  and 
25  is  determined.  This  has  been  changed  to  the  following  steps: 

4 2 

1.  Determine  the  equilibration  of  the  N(  S)  and  N(  D) 
available  at  the  start  of  the  time  interval  using  an 
approximate  form  for  the  electron  density  decay  during 
the  interval.  Nitrogen  atoms  produced  by  the  ion  chem- 
istry are  not  equilibrated  in  this  step. 

? 

2.  Determine  the  loss  of  N(  D)  by  reactions  20  and  22. 

4 

3.  Determine  the  loss  of  N(  S)  by  reactions  19  and  21. 

4.  Determine  the  production  and  loss  of  NO  by  reactions 
5 and  25. 

A o 

5.  Estimate  the  equilibration  of  N(  s)  and  Npp)  after 
the  NO  chemistry  using  the  electron  density  at  the  end 
of  the  time  interval. 
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In  step  1 the  electron  density  is  assumed  to  decay  by  recom- 
bination during  the  time  interval  and  an  effective  recombination  coef- 
ficient determined  from 


- 1 


a = - 
e N 


el 

]Je2 

ciTTTrTiT 


f 4-4) 


where 


N ,,N  - = electron  densities  at  start  and  end  of 
el  c2 

time  interval  respectively. 

4 2 

The  equilibration  of  N(  S)  and  N(  D)  by  reactions  23  and  24  can  be 
determined  by  assuming  that  the  total  number  of  nitrogen  atoms  re- 
mains constant, 

Nt  = [N(4S)]1  ♦ [N(2D) 1 j (4-5) 

and  using  the  assumed  form  for  the  electron  density  decav.  Solving 
2 

for  N(  P)  gives 
2 

AMJ-  . . . k24)  [N(2D)]NC  . k2J  NcNt  (4.6) 


2_  [N C 2D)  ] j - fN(2D)  ]n 


N(“D)2  = 


[1  + ac  Nel(W]y 


[N(^D)1 


(4-7) 


where 


2 2 ^ 

D)L  = k-  "T  V ' NT 
k23  K24  1 


k k 
k23  ♦ 24 

a 


The  concentration  of  N(  S)  after  equilibration  is  then 
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(4-8) 


[N(4S)]2  = nt  - [N(2D)}2  . 

2 

In  step  2 the  loss  of  N(  D)  due  to  reactions  20  and  22  is 
determined.  For  these  reactions 

2 

^011=  - k2Q[N(2D)][02]  - k22[NC2D)  ] [NO]  (4-9) 

= k20[N(2D)][O2]  - k22[N(2D)][NO]  . (4-10) 


These  equations  are  similar  to  the  ozone-oxygen  equations  in  the  D- 
region  for  which  the  following  approximate  solutions  have  been  given 
by  Scheibe  (Reference  4-7)  . 


«"«0 


(4-11) 


n = n0  exp(-A(  1 + £jn)t) 

where 


t = time  interval 
£ = l [NO] 

n = f [n(2d) ] 

A = k20  f°2^ 

B k22 

c - po  - - £n(l  - 50) 

D = 1 + C 


5,  ' «o  * (§  ’ 


L7* 


(4-12) 
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In  the  above  equations  it  is  assumed  that  the  concentration  of  0^  re- 
mains constant.  The  loss  of  can  he  computed  from 


dfO,] 


k20TN(  D) 1 f02l 


(4-13) 


Substituting  for  N(  D)  from  Equation  4-12  and  solving  for  0 gives 


[021  = [02lj  exp 


k20  [N(2n^2  ' W^i 
Ad  ♦ U 


(4-14) 


If  the  value  of  N(  D)  is  larger  than  the  value  of  0 at  the  start  of 

2 ^ 

the  time  interval,  the  loss  of  0 2 and  N(  D)  by  reaction  20  is  computed 
first  from 


d[02l 


k2of°2^(°2)  * 6) 


(4-15) 


where 


6 = (n(2D)j  - [02]j)) 


The  solution  for  Equation  4-15  is 


fO  1 = 

1 2J2  1 - G 


(4-16) 


where 


i°2>>  i | 

[021  * 6exp  |-  k20St|  • 


For  k2Q6t  « 1 


[°2h  ~ 1 k206t 

22  1 + k2of°2ht 


(4-17) 
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(4-18) 


2 

The  amount  of  N(  D)  at  the  end  of  the  interval  is 
fN(2D)l2  = [0212  ♦ 6 

and  the  amount  of  NO  produced  is  just  equal  to  the  loss  of  02 

(NO]  = [02lj  - [02]2  . (4-19) 

2 

Next,  the  loss  of  N(  D)  and  NO  by  reaction  22  is  computed.  The  solu- 
tion is  similar  to  that  given  above  where  the  0.  is  replaced  with  the 
smaller  of  the  N(“P)  and  NO  and  6 is  equal  to  the  difference  between 
them. 

4 

In  step  3 the  loss  of  N(  S)  due  to  reactions  19  and  21  is 
determined  using  the  solutions  described  for  step  2. 

In  step  4 the  production  and  loss  of  NO  by  reactions  5 and 
25  are  completed  using  the  solutions  described  in  Reference  4-2. 

In  steps  2,  3,  and  4 the  appropriate  changes  in  neutral 
species  arc  made  to  maintain  conservation  of  atoms  and  the  energy  re- 
leased by  each  reaction  is  determined  in  order  to  modify  the  pressure 
and  gas  temperature  at  the  end  of  the  time  interval. 

2 

Since  it  is  possible  to  significantly  reduce  the  N(  D)  rela- 

4 

tive  to  the  N(  S)  in  steps  2 and  3,  it  is  necessary  to  estimate  the 

4 2 

equilibration  of  N(  S)  and  N(  D)  to  prevent  oscillation  in  the  pre- 
dicted species  concentrations.  In  performing  this  estimate  the 
electron  density  is  assumed  constant  with  a value  equal  to  that  at 
the  end  of  the  time  interval.  This  prevents  using  electrons  for 
equilibration  that  were  used  in  the  formation  of  the  nitrogen  atoms. 
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SECTION  5 

POST  STABILIZATION  DEBRIS  MODELS 


INTRODUCTION 

In  the  current  WEPH  code  the  debris  cloud  location  and  size 
after  stabilization  are  computed  from  one  of  two  models  selected  as 
a user  input  option.  One  model  neglects  motion  of  the  debris  center 
and  determines  debris  expansion  analytically  in  terms  of  an  empirical 
wind  velocity  profile.  The  second  model  determines  the  location  and 
size  of  the  debris  region  from  a numerical  evaluation  of  debris  mo- 
tion in  atmospheric  wind  fields  specified  as  a function  of  location 
and  time. 

For  debris  stabilization  altitudes  above  about  100  km  the  two 
models  result  in  similar  debris  sizes  and  the  center  of  the  debris 
region  determined  numerically  is  generally  within  the  debris  region 
determined  from  the  analytic  model.  However,  for  debris  stabiliza- 
tion altitudes  below  about  60  km  the  debris  size  determined  numeri- 
cally is  much  smaller  than  that  determined  from  the  analytic  model. 
The  following  is  a brief  review  of  the  two  models  with  suggested 
revisions  to  bring  them  into  better  agreement. 

NUMERICAL  MODEL 

In  the  numerical  model  the  location  of  particles  separated  in 
altitude  between  the  bottom  and  top  of  the  debris  are  followed  in 
the  wind  field  by  computing  their  motion  during  N time  steps  from 

AX  = v AT  (5-1) 

where  AX  is  the  east-west  or  north-south  motion  during  a time  step, 
v is  the  average  east-west  or  north-south  wind  velocity  at  the  start 


of  the  time  step  and  T is  the  duration  of  the  time  step  (tt  1 hour). 

At  the  end  of  N time  steps  the  debris  center  is  determined  as  the 
center  of  the  particles  and  the  debris  radius  is  determined  from  the 
particle  dispersion. 

For  east-west  motion  let 

— T 

X = location  of  debris  center  at  T hours  (AT  = rr ) 

N 

a = standard  deviation  of  debris  center  at  T hours 
x 

If  the  wind  velocities  at  the  start  of  each  time  step  are  uncorrelated, 
then  (see,  for  example  Reference  5-1) 


o 

X 


+ o 


v2 


+ . . . 


where 


(5-2) 


for 


a . = standard  deviation  for  the  velocity  at  the  start  of 
vx  7 

the  ith  time  step 


°vl  = av2  - °VN  " 0v 
a T 

a = ^ a AT  * 
x v 

If  the  wind  velocities  at  the  start  of  each  time  step  are  correlated, 
then 

°x  = + av2  + “*  avN^T  (5_3) 

and  for  equal  standard  deviations 

a = No  AT  = a T . (5-4) 

X V V 

In  Reference  5-2  variations  of  wind  velocities  from  the  mean 
are  described  in  terms  of  small  scale  and  large  scale  components.  The 
large  scale  component  is  related  to  tides,  gravity  waves  and  planetary 
waves,  and  the  small  scale  component  is  related  to  turbulence.  In  the 


- 


V* 


horizontal  direction  the  small  scale  correlation  distances  vary  from 
about  20  km  at  the  surface  to  about  150  km  at  100  km  altitude  and  the 
large  scale  correlation  distances  vary  from  about  900  km  at  the  sur- 
face to  about  1500  km  at  100  km  altitude.  The  small  and  large  scale 
standard  deviations  are  related  to  the  total  standard  deviation  by 
(see  Reference  5-2) 

2 2 
< - fL  °T 


where 


(i  - f l)<4 


= total  standard  deviation 


a^Oj,  = large  and  small  scale  standard  deviations 
respectively. 

The  fraction  f^  of  the  total  variance  contained  in  the  large  scale 
variance  is  given  in  Reference  5-2  as  a function  of  altitude  and  lati- 
tude and  varies  from  about  0.85  at  20  km  altitude  to  0.65  at  100  km 
altitude. 

Assuming  that  the  step  sizes  used  in  computing  the  debris  lo- 
cation are  large  in  comparison  to  the  small  scale  correlation  distance 
and  small  in  comparison  to  the  large  scale  correlation  distance,  o 
becomes 


o = i / 0. 25(a  \ + a \ + 


, , + a _ 
vl  v2 


/ °-75Cavl  + °v2  + 


a AT  + 
vN 


a .,)  AT 
vN 


(5-5) 


where  f,  has  been  chosen  as  0.75  and  a . is  the  total  standard  devia- 
L vi 

tion  at  the  start  of  the  ith  step.  For  equal  standard  deviations 

(5-6) 


V 

a = 0.5  — ♦ 0.87  a T 

x /F  v 


99 


Thus,  since  there  is  a high  degree  of  correlation  in  the  wind  velocity 
in  regions  of  about  1000  km,  the  standard  deviation  of  the  *ebris  cen- 
ter is  essentially  the  same  as  for  perfect  correlation.  This  results 
in  a considerably  larger  uncertainty  in  determining  the  location  of 
the  debris  centers  than  indicated  in  Reference  5-1  where  uncorrelated 
velocities  were  assumed. 

The  debris  radius  in  the  numerical  model  is  found  in  terms  of 
the  particle  displacement.  The  east-west  debris  diameter  is  found 


Df:W  “ ' Xa 


(5-7) 


where 


= location  of  the  two  particles  having  the  greatest 
east-west  separation  after  N time  steps. 

While  this  is  a reasonable  estimate  of  the  debris  size  when  the  hori- 
zontal velocities  at  different  altitudes  are  correlated,  it  can  under- 
estimate the  size  when  the  velocities  are  uncorrelated.  In  Reference 
5-2  the  small  scale  vertical  correlation  distance  is  given  as  about 
4 km  and  the  large  scale  correlation  distance  as  about  30  km  (mid- 
latitudes).  Thus,  for  typical  debris  regions  the  particles  are  sepa- 
rated by  a distance  comparable  to  or  larger  than  the  small  scale 
correlation  distance  and  less  than  the  large  scale  correlation  distance. 


‘ Xa  • *b 


Xa,Xb  = mean  position  of  Xg  and  X^ 

o ,o  , = standard  deviations  of  the  uncorrelated  variation 

Xa  XD 

in  Xg  and  X^ 


Z = Xa  - *b 


Xa  > *b 


(5-9) 


100 


r 


i 


rr ~ r 

a,  = /a  + a , 

Z xa  xb 


(5-10) 


Assuming  a normal  distribution  for  Z the  mean  value  of  Dew  is 
derived  as  (see  Appendix  C) 


EW 


where 


F(X) 


Z 1 - 2F 


■*/ 


( 4)1' 


(5-11) 


dt 


Note  that  for  Z 

2a„ 


D, 


™ ' /SF  * 


(5-12) 


In  order  to  estimate  the  mean  debris  spread  from  velocity 
data  (standard  deviation  data  are  not  currently  modeled  in  the  WEPH 
code)  the  variance  in  velocity  is  assumed  to  be  the  same  at  each  alti- 
tude within  the  debris  region  and  equal  to  the  value  of  the  largest 
velocity. 


a . - v , 
vi  mi 


(5-13) 


where  vmi  is  the  largest  velocity  at  the  start  of  the  ith  step.  Then, 


a_  = a 
Z x 


°x  ■ °-25(vJl  vmN>'4T>2 

where  as  before  the  quantity  1 - f^  is  taken  as  0.25. 
east-west  diameter  becomes 


(5-14) 

(5-15) 

The  average 
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m 


I 


t 


I 

l 


I 

'5 


dfw  a z l - 2F 


(- 1 )) 


Tf/I"* 


2 

Vm2  + 


V M ^Te 

mN 


■m 


(5-16) 


or  considering  the  approximations 

* z + 4f  / 2 + 2 


EW 


ml 


v - + ...  v ..  ATe 
m2  mN 


(5-17) 


A similar  calculation  is  made  for  the  north-south  debris  di- 
ameter and  the  debris  radius  computed  from 


R - R + ifeu 

R Rs  2 Tj 


(5-18) 


where 


R 


S = debris  radius  at  stabilization  time. 


Figure  5-1  shows  comparisons  of  the  debris  radii  computed  from  the 
current  model  and  from  the  above  relations. 

ANALYTIC  MODEL 

In  the  analytic  model  the  debris  radius  is  computed  from 

(5-19) 


R - Rs + VT  - V 


where 


v min 
wp 


[■T-] 


wp 


= characteristic  wind  speed  at  the  debris 
stabilization  altitude 


DV(TS)  = vertical  dimension  of  the  debris  at  the 
stabilization  time 
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T * Tj  (hours) 


Figure  5-1.  Comparison  of  current  and  revised  numerical  models 


H = the  characteristic  distance  over  which  the 

wp 

wind  either  changes  direction  or  velocity 
sufficiently  to  disperse  the  debris,  taken 
as  10  km. 

The  profile  of  the  wind  speed  is  assumed  to  be 


0.01 

+ 0.02 

"D 

15 

km/sec 

0 < hp  < 15 

0.05 

- 0.02 

hD- 

15 

15  < hR  < 30 

nr 

0.01 

♦ 0.04 

HP 

30 

30  < hp  < 80 

5d 

o.os 

hD>8° 

(5- 


where  h„  is  the  debris  stabilization  altitude. 


The  above  model  predicts  much  larger  debris  regions  at  low 
altitudes  than  the  numerical  model  due  to  overestimating  the  amount 
of  wind  shear  and  neglecting  the  difference  in  east-west  and  north- 
south  velocities.  In  oi'der  to  reduce  the  debris  expansion  the  follow- 
ing revised  model  is  proposed.  For  times  after  debris  stabilization 
compute  the  debris  radius  from 

(5-2” 


where 


vw  = effective  east-west  debris  spreading  velocity 
rw  = ratio  of  the  north- south  to  east-west  wind  velocity. 
The  proposed  debris  spreading  velocity  is 


where 


v = v min 
w wp 


, r wi 

T %r\ 


C5-22) 


wp 


0.01  km/sec 

- 0.11  ♦ 0.002  hr 


n 


0.05 


hn  s 60 

60  < hn  < 80 

hD  >80 


and  Dy(Tg)  and  are  as  previously  defined.  The  proposed  velocity 
ratio  is 


rw  = 


hp<6° 

60  < hp  < 80 
hp  > 80  . 


(5-25 


Figure  5-2  shows  comparisons  of  debris  radii  computed  from  the  current 
and  revised  analytical  models.  Also  shown  on  Figure  5-2  is  the  debris 
radii  for  the  revised  numerical  model.  The  proposed  changes  result  in 
nearly  equal  debris  expansions  for  the  analytical  and  numerical  models. 
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APPENDIX  A 

REACTIONS  AND  REACTION  RATE  COEFFICIENTS 
USED  IN  D-REGION  CHEMISTRY  MODEL 


The  following  reactions  and  reaction  rate  coefficients  are 
for  use  in  the  D-region  chemistry  model  (see  Section  4).  The  reac 
tion  rate  coefficients  given  arc  in  the  form 


Reactions 


Neutral  Species  Reactions: 

1 Oj  ♦ hv  ♦ ♦ Oj  ♦ 0 

2 03  ♦ hv  - 02(’a)  ♦ 0 

3 02  ♦ hv  -*0*0 

« N0?  ♦ hv  - NO  ♦ 0 

5 N(z0)  ♦ 0,  - NO  ♦ 0 


6 0^1  A)  Op  ♦ hv 

7 02('a)  ♦ N - 02  ♦ N 

8 02(’a)  ♦ N2-  0z  ♦ N2 

9 o2(1a)  ♦ 0 - o2  ♦ 0 

10  0,(’a)  ♦ 0,-  20, 


Current  Code  Values 


(Effective  values  computed 
uithin  model) 


Revised  Code  Values 


7.0[-l2]  .0.5 


16  N ♦ 0 

17  N ♦ 0 ♦ H 

18  N ♦ 02 

19  N ♦ 03 

20  N ♦ 0 


- NO 

- NO  ♦ M 
-*  NO  ♦ 0 

- NO  ♦ 02 
■*  NO.  ♦ hv 


1.21-17] 
1 - 1 1*  32] 
3.3[-12] 
3.«[-U] 
6.6[-17] 


l.«t-17] 

U 

U 

3.1[-11] 

u 


Read  7.0[-12]  as  7.0  * 10"2 
Rate  coefficient  unchanged  from  current  code  value 


Reactions 


Current  Code  Values  Revised  Code  Values 

A B C A B C 


41 

NOj  ♦ H02 

- HN02  ♦ 02 

3. 0 [- 1 4] 

42 

OH  ♦ H0? 

- h?o  ♦ o2 

2.0[-11] 

43 

o('a)  ♦ h2o  - oh  a oh 

2.3[-10] 

44 

Of1.*.)  ♦ N2 

- 0 ♦ n2 

2.5[-ll] 

45 

0( 'a)  a 02 

- 0 ♦ o2 

3.0[-11] 

40 

o(’a)  ♦ H2  - h a oh 

1.0t-10] 

4; 

0(’i)  ♦ H202  ♦ oh  ♦ ho2 

5.2[-10] 

48 

H ♦ H02 

• H2  * °2 

4.2[-ll] 

350 

49 

H ♦ H02 

- h?o  ♦ o 

8.3[-11] 

500 

50 

H ♦ H02  - OH  a OH 

4.2[-10] 

950 

Positive  Ion  and  Electron  Reactions: 


100 

I 

HjO*  ♦ H20  ♦ H 

•• 

H*(H20)2  ♦ M 

3.4[-77] 

-2.0 

u 

1 

H*(H20)2  a H 

• 

H30*  ♦ H20  ♦ N 

8.0 

-2.0 

18000 

4.01-1] 

-5.0 

17100 

102 

HjO*  ♦ N2  a M 

-• 

H30*(N2)  V M 

1.4C-30] 

-2.0 

1 . 4[- 30] 

-4.0 

103 

HjO*(N2)  a M 

- 

HjO*  ♦ N2  ♦ M 

1.2[-8] 

-2.0 

2780 

2.0[-6] 

-5.0 

4500 

104 

h3o*(ho)  ♦ h2o 

- 

H*(H20)2  a HO 

1.4[-9] 

U 

105 

HjO*(HO)  ♦ N 

- 

HjO*  ♦ HO  ♦ * 

l.OC-1] 

-2.0 

11800 

5.0[1] 

•5.0 

13000 

10SA 

h3o*(n2)  ♦ h2o 

- 

H*(H20)2  ♦ N2 

1.0[-9] 

u 

106 

H*(H20)2  ♦ H20  ♦ M 

H*(H20)3  ♦ M 

2. 3[-27] 

o 

N 

1 

2.3C-27] 

-4.0 

107 

H*(H20)j  ♦ M 

H*(H20)2  ♦ H20  ♦ N 

i.oc-i] 

© 

N 

l 

11200 

8.0C-2) 

-5.0 

11000 

108 

! 

H+(H20)j  ♦ H20  a M 

■* 

H*(H20)4  a M 

2.  *[-27] 

-2.0 

2.4C-27] 

-4.0 

109 

Hf(H20)4  ♦ M 

■* 

H*(H20)j  a H?0  * 4 

i.oc-i] 

-2.0 

8600 

5.6C-2) 

-5.0 

8360 

110 

h*(h2o)4  ♦ h2o  ♦ m 

h*(h2o)5  ♦ n 

8.8[-28] 

-2.0 

9.0[-28] 

-4.0 

— 

H*(H20)5  ♦ H 

♦ 

H*(H20)4  a H?0  a H 

3.0[-3] 

-2.0 

6000 

2.81-1] 

-5.0 

7700 

112 

N*  • 02 

Oj  a N(*0) 

3.0C-10] 

U 

113 

N*  a 02 

-♦ 

NO*  ♦ 0 

3.0C-10] 

U 

114 

NO*  ♦ C02  ♦ »4 

N0*(C02)  a N 

2.0C-29] 

-2.0 

7.0[-»] 

-3.0 

115 

N0*(C02)  ♦ M 

NO*  a C02  a H 

1.0C-8] 

•2.0 

5000 

3.8C-6] 

-4.0 

4590 

no 

_ ’ j 


116 

NO*  ♦ H2  0 ♦ M 

N0*(H2O)  ♦ M 

1.51-28] 

-2.0 

1.8[-28] 

-4.7 

1 

117 

NO*  ♦ N?  ♦ M 

NO*(N2)  ♦ M 

2.0[-31] 

-2.0 

2. 0[- 31 ] 

-4.4 

118 

no*(n2)  ♦ m 

NO*  ♦ N2  ♦ M 

1.2[-B] 

-2.0 

2780 

6.3[-8] 

-5.4 

2450 

'>  119 

no*(co2)  ♦ H20 

N0*(H20)  ♦ C02 

1.0[-9] 

U 

120 

N0*(H20)  ♦ HO 

h3o*  ♦ no2 

1.01-9] 

6.0[-ll] 

1 ~ 

N0*(M20)  ♦ H02 

HjO*  ♦ NO  ♦ 02 

1.0t-9] 

l.Ot-10] 

122 

N0*(H2O)  ♦ H?0  ♦ M * 

N0*(H20)2  ♦ M 

l.lt-27] 

•2.0 

l.H-27] 

-4.7 

123 

N0*(H20)2  ♦ h 

N0*(H20)  ♦ H20  ♦ M 

1.01-4] 

-2.0 

7000 

9. 7[-3] 

-5.7 

8100 

124 

N0*(H20)2  ♦ h2o  ♦ H - 

N0*(H20)3  ♦ M 

1.6[-27] 

-2.0 

1.0[-27] 

-4.7 

125 

N0*(H20)j  ♦ m 

NO*(H20)2  ♦ H20  ♦ M 

2.0[-2] 

-2.0 

7000 

1.2t-2] 

-5.7^ 

•^6800 

126 

N0*(H20)3  ♦ h2o 

- 

h*(h20)3  ♦ hno2 

7.0[-11] 

U 

/ 

y 

127 

NO*(N2)  ♦ C02 

-* 

N0*(C02)  ♦ n2 

1.01-9] 

u 

128 

NO*(N2)  ♦ HjjO 

■* 

N0*(H20)  ♦ N2 

1.0[-9] 

u 

129 

NO;  ♦ NO 

-* 

NO*  ♦ N02 

2.9[-10] 

u 

130 

N*  ♦ 0, 

-♦ 

°2*N2 

S.0[-11] 

-0.5 

5.0[-11] 

-0.8 

131 

0*  ♦ N2 

•* 

NO*  ♦ N 

1 . 2[- 12] 

-0.5 

1.2C-12] 

-1.0 

132 

0*  ♦ N2  ♦ M 

- 

NO*  ♦ N ♦ M 

!.6[-29] 

-2.0 

e.oc-29] 

•2.0 

133 

o*  ♦ o2 

■* 

oj  ♦ 0 

2-0[-11] 

-0.5 

2 -0[- 11 ] 

-0.4 

134 

0*  ♦ HjjO  ♦ M 

-* 

0*(H20)  ♦ M 

2.8[-28] 

-2.0 

U 

135 

oj  ♦ N 

NO*  ♦ 0 

1 . 8[- 10] 

1 * 2[- 10] 

Ill 


Reactions 

Current  Code  Values 
AS  C 

Revised  Code  Values 
ABC 

136 

0^  ♦ NO 

NO*  ♦ 02 

6.3[-10] 

4.5[-10] 

137 

0*  ♦ N0? 

NoJ  ♦ 02 

6.6[-10] 

U 

138 

02  ♦ no2 

- 

NO*  ♦ Oj 

101-11] 

U 

139 

°2  * *2 

NO*  ♦ NO 

1.0t-16] 

u 

140 

oj  ♦ N2  ♦ N 

■* 

02{N2)  ♦ N 

8 . 0[ - 3 1 ] 

-2.0 

9.0[-31] 

-2.0 

141 

02(N2)  ♦ M 

* 

02  ♦ N2  ♦ N 

2 . Of  - 1 1 ] 

U 

142 

oj  ♦ 0?  ♦ N 

■* 

Oj  ♦ N 

2.8[-30] 

-2.0 

3.9[-30] 

-3.2 

143 

oj  ♦ M 

- 

Oj  ♦ 02  ♦ N 

5.0[-7] 

-2.0 

5000 

2.0[-5] 

-4.2  5400 

144 

02(H20)  ♦ HjO 

* 

HjO  ♦ HO  ♦ 02 

2.0[-10] 

U 

115 

o*(h2o)  ♦ h2o 

HjO*{HO)  ♦ 02 

1.0[-9] 

U 

146 

02(H20)  * NO 

■* 

NO*  ♦ H20  ♦ 02 

1.0[-10] 

U 

147 

02(h20)  ♦ 02(’a) 

-♦ 

o2  ♦ h2o  ♦ oz 

1.0t-10] 

U 

148 

oJ(n2)  ♦ 02 

«* 

1.0[-9] 

U 

149 

oj  ♦ H20 

- 

0j(Hz0)  ♦ 02 

)5[-9] 

U 

150 

ojoyi)  ♦ o2 

-* 

oj  ♦ H20 

2.0[-10] 

2300 

U 

151 

oj  t NO 

-* 

NO*  ♦ 202 

S.0C-10] 

U 

152 

oj  ♦ N02 

NO*  ♦ 20? 

5.0[-10] 

U 

153 

oj  t 0 

- 

°5*°3 

3.01-10] 

u 

154 

h4(h20)j  ♦ e 

-* 

H ♦ 3H20 

4.0[-6] 

5.1  [-6] 

155 

h+(h20)4  ♦ e 

H t 4HjO 

4.9[-6] 

6.11-6] 
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Ructions 


156 

H*(M20)5  ♦ e 

■* 

H ♦ 5H?0 

157 

NO*  ♦ e 

-♦ 

J«(2o)  ♦ \h 

158 

N0*(H20)2  ♦ e 

♦ 

2H?0  ♦ NO 

159 

N0*(H20)3  ♦ e 

-* 

»20  ♦ NO 

160 

N0*(N2)  ♦ e 

-♦ 

NO  ♦ N? 

161 

n!  ♦ e 

-♦ 

N(*t»  ♦ N 

162 

02  ♦ e 

♦ 

o{'o)  ♦ 0 

163 

02{H20)  ♦ e 

-* 

«2o  ♦ o2 

164 

oj  ♦ e 

♦ 

202 

165 

e ♦ NO, 

♦ 

no: 

I 

I 

Current  Code  Velues 

devised  Code  Velues 

A 

« 

C 

A 

9 C 

6.0[-6] 

7.4C-6] 

«.0[-7] 

-1.0 

U 

2.0[-6] 

-1.0 

2.0[-6] 

-0.2 

3.0[-6] 

-1.0 

5.0[-6] 

-0.2 

1.0[-6] 

-1.0 

1-51-6] 

2. 7[-7] 

•0.2 

2.71-2] 

2-U-7] 

-0.7 

2.U-7] 

•0.63 

1.5[-6] 

•1.0 

1.5C-6] 

-0.2 

i 

2.0[-6] 

-1.0 

U 

I 

4.01-11] 

U 

I 

1.0C-31J 

U 

1*1-12] 

1.5 

4990 

u 

1.4{-29] 

-1.0 

600 

u 

2.7C-10] 

0.5 

5590 

1.9[-10] 

0.5  5600 

*.0[-12] 

1.5 

U 

2.2C-T0]  U 


2.01-10]  2.51-10]  -o.a 

2.0[-10]  U 

J.0t-I0]  U 

J.0[-10]  u 
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I 


1 


' 

I 


Deletions 

Current  Code  Velues 

A B C 

Revised  Code  Velues 

A B C 

I 

j 

176 

0*  ♦ 0 

♦ 

e ♦ Oj 

3.0[-10] 

1 • 5[- 10  ] 

1 

• 

\ 

177 

°j>  ♦ °2<’v 

♦ 

e ♦ 202 

2.0[-10] 

U 

f 

178 

0 ♦ e 

-* 

0 ♦ hv 

1.31-15] 

u 

1 

179 

« ♦ 02  HjO 

* 

Oj  ♦ H20 

1.4[-29] 

u 

I 

Negetivc  Ion  Deletions: 

200 

COj  ♦ H20  ♦ H 

-* 

CO j(  HjO)  ♦ M 

1.0[-28] 

-1.0 

u 

201 

COj  ♦ NO 

NOj  ♦ C02 

1.8[-U] 

1.11-11]  -1.1 

! 

202 

COj  ♦ no2 

* 

NOj  ♦ C02 

8.0t-11] 

2.0[-10] 

203 

COj  ♦ 0 

-* 

oj  ♦ co2 

8 . 0[ - 11] 

1.11-10] 

204 

C0j[H20)  ♦ NO 

- 

NO j<  H?0 ) * C02 

1. •[•!!] 

7. 0 [- 1 2] 

! 

205 

COj  ♦ H20  ♦ M 

- 

coj(H2o)  ♦ n 

5.0[-29] 

•1.0 

u 

206 

COj  ♦ NO 

- 

NOj*  ♦ C02 

4.8[-11] 

u 

i* 

j 

207 

COj  ♦ 0 

-* 

COj  ♦ 02 

i.st-io] 

u 

208 

COj  ♦ Oj 

Oj  ♦ C02  ♦ 02 

1.3C-10] 

u 

209 

COj(HjO)  t NO 

* 

NO j( h2° ) ♦ C02 

5.0[-10] 

1.0[-11] 

210 

NOj  ♦ HjO  ♦ M 

* 

N0j<H20)  ♦ M 

1 . 3[*28] 

-1.0 

u 

i 

NOj  ♦ N02 

- 

NOj  ♦ NO 

4.0C-12] 

2.0[-13] 

212 

NOj  ♦ Oj 

NOj  ♦ 02 

1.8C-11] 

9.0[-ll] 

213 

NOj(H20)  ♦ Oj 

* 

NOjtHjO)  ♦ Oj 

5.0C-10] 

1.01-11] 

214 

NOj  ♦ HjO  ♦ M 

N0j(H20)  ♦ N 

7.5[-29] 

-1.0 

u 

215 

0*  ♦ C02  ♦ H 

•* 

COj  ♦ N 

3.U-28] 

-1.0 

u 

s 


■ 

, 
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1 


Ructions 

Current  Code  Values 

A B C 

Revised  Code  Values 

A B C 

216 

O' 

♦ no2 

♦ 

NOj  ♦ 0 

1.2[-9] 

0 

217 

0* 

♦ 0?  ♦ M 

♦ 

Oj  ♦ M 

l.lf-30] 

-1.0 

« 

218 

0" 

♦°J 

* 

Oj  ♦ 0 

5.3[-10] 

U 

219 

°2 

♦ C02  ♦ N 

■* 

coj  ♦ M 

2.0[-29] 

-1.0 

220 

♦ H20  ♦ M 

•* 

02(H20)  ♦ M 

3.0[-28] 

-1.0 

u 

1 

221 

Oj  ♦ N 

■* 

O'*  NO 

1.0[-10] 

U 

222 

o*  ♦ no2 

* 

NOj  ♦ 02 

a.ot-io] 

1.2C-9] 

223 

oj  ♦ 0 

-* 

o'*  o2 

3. 3[- 10] 

1.5[-10] 

224 

Oj  ♦ 02  ♦ M 

• 

Oj  ♦ M 

3.5[-31]  -1.0 

U 

225 

Oj  ♦ M 

Oj  ♦ 02  ♦ M 

1.0[-3]  -1.0 

7500  2.0[-5]  -1.0  6300 

226 

°2  * °3 

* 

oj  * o2 

4.0C-10] 

6. Of- 1 0] 

227 

0*(HZ0)  ♦ C02 

♦ 

coj  ♦ h2o 

5.8C-10] 

U 

228 

02(H20)  ♦ NO 

<* 

NOj*  ♦ H20 

3.1C-10] 

u 

229 

0'(H20)  ♦ Oj 

♦ 

oj  ♦ m20  ♦ 02 

3.1[-10] 

2.3[-10] 

230 

Oj  ♦ C02 

* 

coj  ♦ 02 

5.5C-10] 

U 

231 

Oj  ♦ NO 

* 

NOj  ♦ 02 

i.oc-11] 

2. 8 [- 1 2] 

232 

oj  ♦ no2 

* 

NOj  * 0, 

1.41-10] 

0 

233 

Oj  ♦ N02 

-• 

NOj  ♦ 02 

1.4[-10] 

2.8C-10] 

234 

Oj  ♦ 0 

-* 

°j*02 

i.ot-ii] 

3.2[-10] 

235 

Oj  ♦ C02 

coj  ♦ 02 

4.3C-10] 

U 

115 


| •*«(*•>*  f 
• I-?'-**! 

IS 


Ructions 


Current  Code  Values 
ABC 


Revised  Code  Values 
A B C 


236 

CO*  ♦ o2 

♦ 

o;  ♦ co2 

4. 3[-l0] 

3000 

8.41-12] 

237 

Oj  ♦ H?0 

<* 

0*(H20)  ♦ 02 

U 

238 

Oj  ♦ NO 

- 

NO**  ♦ 02 

2.5[-10] 

0 

239 

o:  ♦ o 

♦ 

°3  * °2 

4.01-30] 

0 

240 

°4  * °3 

* 

0*  ♦ 202 

5.0[-10] 

3.0[-10] 

241 

NO**  ♦ NO 

■* 

NO*  ♦ N02 

1.5C-11] 

u 

242 

H*  ♦ H* 

* 

Products 

2.0[-7]  -0.5 

6. 0[-8] 

243 

H*  ♦ H*  ♦ M 

♦ 

Products 

5.6[-25]  -1.5 

1.0[2S] 

! 


APPENDIX  B 


REACTIONS  AND  REACTION  RATE  COEFFICIENTS 
USED  IN  E-  AND  F-REGION  AND  HEATED 
REGION  CHEMISTRY  MODELS 


The  following  reactions  and  reaction  rate  coefficients  are 
used  in  the  E-  and  F-region  chemistry  model  (see  Section  4)  and  the 
heated  region  chemistry  models  (see  Reference  4-1).  The  reaction 
rate  coefficients  given  are  in  the  form 

-c/T 


k ■ k{m)  e 


Ructions 

Current  Code  Values 

A B C 

Revised  Code 
A B 

Values 

C 

1 

n*  ♦ o2 

♦ 

NO*  ♦ 0 

3-Of-IO]* 

If 

2 

n*  ♦ o2 

-♦ 

02  ♦ N 

3.01-10] 

U 

3 

n*  ♦ o2 

♦ 

0*  ♦ NO 

1 .0[-12] 

U 

4 

N*  ♦ N2 

N2  ♦ N 

4.0[-11] 

20060 

4.0[-U] 

20100 

5 

N*  ♦ NO 

♦ 

NO*  ♦ N 

8.0[-10] 

U 

6 

N*  ♦ 0 

•» 

0*  ♦ N 

1.0[-12] 

2.6t-13] 

1.0 

7 

N*  ♦ E 

N ♦ hv 

4.4C-12] 

-0.75 

U 

8 

N*  ♦ E ♦ E 

♦ 

N ♦ E 

1.2H9] 

-5.0 

U 

9 

o*  ♦ o2 

- 

0*  ♦ 0 

2.0[-11] 

-0.5 

2.0C-11] 

10 

o*  ♦ n2 

NO*  ♦ N 

6.0[-14] 

2.0 

-900 

u 

11 

0*  ♦ E 

0 ♦ hv 

<•<[-12] 

-0.75 

u 

12 

0*  ♦ E ♦ E 

* 

0 ♦ E 

1.2[-19] 

-5.0 

u 

13 

NO*  ♦ E 

♦ 

N(4S)  ♦ 0 

i.oc-73 

-1.0 

u 

14 

NO*  ♦ E 

N(2D)  ♦ 0 

3.0C-73 

-1.0 

0 

IS 

N(4S)  ♦ 0 

NO*  ♦ E 

2.1[-13] 

0.5 

31900 

2.  It*  13] 

0.5 

31900 

N ♦ N - N?  ♦ hv 
N(4$)  ♦ 02  * NO  ♦ 0 
N(2D)  ♦ 0, 


♦ NO  ♦ 0 


Reactions 


Current  Code  Values  Revised  Code  Values 

A G C ABC 


21 

N(4S)  ♦ NO 

- N2  ♦ 0 

2.71-11] 

U 

22 

N(ZD)  ♦ NO 

♦ N2  ♦ 0 

7.01-11] 

6.01-11] 

0.5 

23 

N(4S)  ♦ E 

- N(20)  ♦ E 

1.3C-9] 

26000 

1. 3[-9] 

!7700 

24 

N(2D)  ♦ £ 

- N(4$)  ♦ E 

5.0C-10] 

u 

25 

o ♦ h2 

- NO  ♦ N 

2.5[-10] 

37900 

1.21-10] 

37800 

26 

E ♦ 0?  ♦ 02 

* °2  * °2 

1.4[-29] 

-1.0 

600 

u 

27 

°2  * °2 

- E ♦ 02  ♦ 02 

2.7[-10] 

0.5 

5590 

1.91-10) 

0.5 

5600 

28 

O’  ♦ 0 

•*  Oj  ♦ E 

3.0[-10] 

1.5[-10] 

29 

°3  * E 

- O’  ♦ 0 

1.41-3] 

-1.98 

7194 

1.5[-3] 

-2.5 

7200 

» 

0'  ♦ 0?  ♦ N 

♦ Oj  ♦ N 

3.5C-31] 

-1.5 

3. 5[- 31  ] 

-1.0 

31 

oj  ♦ n 

- Oj  ♦ 02  ♦ N 

2.U-4] 

-1.7! 

6230 

2 . 0[-5] 

-1.0 

6300 

32 

0*  ♦ 0 

- oj  ♦ o2 

4.0C-10] 

0 

33 

°3  * °2 

« Oj  ♦ 0 

4.4(-U] 

-0.23 

25423 

5. 7[-10] 

-2.0 

24200 

34 

O’  ♦ NO 

- OONO*  ♦ 02 

2.5t-10j 

0 

35 

OONO’  ♦ 0? 

- Oj  ♦ NO 

1.8[-8] 

-0.97 

28338 

2.5[-9] 

-0.25 

22900 

36 

0’  ♦ 0 

* °2  * °2 

1.0[-11] 

2.5C-10] 

37 

°2  * °2 

- Oj  ♦ 0 

5.0[-ll] 

0.15 

27711 

5.0[- 12  ] 

29000 

38 

O’  ♦ NO 

■»  NOj  ♦ 02 

1.01-11] 

U 

39 

02  ♦ NOj 

- Oj  ♦ NO 

6.2C-U] 

-0.09 

28281 

5.5[-ll] 

28100 

40 

OONO*  ♦ NO 

* NOj  ♦ NOj 

1.51-11] 

U 
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Reactions 

Current  Code  Values 

A B C 

Revised  Code  Values 
ABC 

61 

AlO*  ♦ M 

- 

At02  ♦ E ♦ M 

7. 3[-9]  -0.03 

51429 

1.01-9] 

50800 

62 

AlO*  ♦ £ 

- 

At  ♦ 0 

5.0[-7]  -1.0 

U 

63 

0 ♦ N0Z 

- 

NO  ♦ 02 

9.0[-12] 

u 

64 

NO  ♦ 02 

- 

o ♦ no2 

6.4[-12]  -0.25 

23200 

1.5[-17] 

-0.25 

23200 

6S 

02  ♦ NO 

- 

NO"  ♦ 0 

1.01-12] 

U 

66 

no2  ♦ 0 

0'  ♦ NO 

6.2[-6]  -1.5 

3900 

1-2C-9] 

-1.5 

3900 

67 

0 ♦ NO  ♦ M 

-♦ 

N02  ♦ M 

Not  Included  In 
Current  Model 

1-61-32] 

-580 

68 

NO  ♦ Oj 

- 

N02  ♦ 02 

2.11-12] 

1450 

69 

0 ♦ N? 

- 

NO  ♦ N(4S) 

1.23[-10] 

37800 

70 

o ♦ h2o 

- 

OH  ♦ OH 

9. 5 [- 1 1 ] 

9000 

71 

n(4s)  ♦ no2 

- 

N?0  ♦ 0 

2.0[-ll] 

800 

72 

n(4s)  ♦ o3 

NO  ♦ 02 

3.11-11] 

1200 

73 

0 ♦ 02  ♦ M 

* 

03  ♦ M 

l.U-34] 

-510 

74 

N2  *°2 

°2  * N2 

5. 0[- 1 1 ] 

-0.8 

75 

N2  ♦ NO 

- 

NO*  ♦ N2 

3.3C-10] 

76 

N(4S)  ♦ 0 

- N(Z0)  ♦ 0 

5.0E-13] 

0.50  27700 

77 

N(ZD)  ♦ 0 

- N(4S)  ♦ 0 

2.0[-13] 

0.50 

78 

0 ♦ N2  ♦ H 

- N20  ♦ M 

1.0[-34] 

7500 

79 

o ♦ n2o 

- NO  ♦ NO 

17[-10] 

14000 

80 

0 ♦ OH 

* H ♦ 02 

4.2[-11 ] 

Reactions 

Current  Code  Values 
ABC 

Revised 

A 

Code  Values 

B C 

81 

H ♦ 0? 

• OH  ♦ 0 

Not  Included  In  Current 
Model 

l.l[-9] 

■0.30 

8400 

82 

N ♦ NO  ♦ M 

♦ NOj  ♦ M 

l.Ot-31] 

2.50 

83 

h*o3 

• OH  ♦ 02 

2.6[.U] 

84 

Oj  ♦ M 

» 0 ♦ 02  ♦ M 

3.8[-8] 

1.50 

12200 

85 

0*  ♦ NO 

- NO*  ♦ 02 

6.3C-I0] 
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APPENDIX  C 

PROBABILITY  DISTRIBUTION  OF  THE  ABSOLUTE 
DIFFERENCE  OF  TWO  INDEPENDENT  RANDOM  VARIABLES’ 


Let  Xj  and  X2  be  two  independent  random  variables  with  mean 
values  of  and  ar>d  standard  deviations  of  and  c^.  Define  a 
new  variable  y by 


y = |x.  - x 


1 2 1 


The  following  deviation  is  for  the  mean  value  and  standard  deviation 
of  y.  First  define  z = x^  - x^.  Then  for  > u^ 

yz  = U1  ‘ ^2 

f~2  2 

°z  = \/°  1 + °2 


If  we  assume  that  z is  normally  distributed 


P[z(t)J 


- 

1 2 1 O ) 

2 ; 


then  y(t)  = |z(t)|  has  the  distribution 
P[y(t)]  = P[z(t)]  + P[z(-t)] 


v/i7o. 


iWL-K 


-t  - U 


The  mean  value  of  y is  found  from 


t > 0 


Prepared  by  James  H.  Thompson 
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Let  s 
then 


124 


] 

The  cumulative  distribution  function  of  a standardized  normal  random 
variable  is 

, r -h2 

F(x)  = -4=  / e dt 

V 2tt  J 

V -00 

The  values  of  F(x)  are  listed  in  tables  or  can  be  found  from  the  error 
function  using  the  relation 


I.-. 

t r ■ 


JSL 


\ 


This  is  in  the  standard  form  for  a normal  random  variable  so  that 

~ 2 2 2 , 2222 

t = oz  + yz  , hence  °y  = °z  + " yy 

Hence  for  the  variable  y = | = x2l 


where 


y = y,  - y»  and  we  have  assumed  y.  > y- 
2 12  i — 2 


j 

I 


I 
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